Explicit Construction of the 
Representation of the Braid Generator a 
Associated with the One-Parameter Family of 
Minimal Typical Highest Weight (0, | a) 
Representations of Uq[gl{2\l) 
and its use in the Evaluation of the 
Links-Gould 
Two- Variable Laurent Polynomial Invariant 
of Oriented (1,1) Tangles 

David DeWit 
BSc, BScAppGeophHons 
BSc, PGDipSci, MScSt, BScHons 



A thesis submitted in the 
Department of Mathematics 
at 

The University of Queensland 
for the requirements of the degree of 
Doctor of Philosophy 

25 November 1998 



Statement of Originality 



The work presented in this thesis is, to the best of my knowledge and behef, original 
and my own work, except as acknowledged in the text, and the material has not been 
submitted, either in whole or in part, for a degree at this or any other university. 

The material presented herein will be published in: 



David De Wit, Louis H Kauffman and Jon R Links. An Invariant of Knots 
and Links Brought to You by Uq[gl{2\l)]. To appear in the Proceedings of 
the XXII International Colloquium on Group Theoretical Methods in Physics 
(Group 22), held at the University of Tasmania, Hobart from 13-17 July 1998. 
International Press. 



[ p.3|| David De Wit, Louis H Kauffman and Jon R Links. On the Links-Gould 
Invariant of Links. Accepted for publication by the Journal of Knot Theory 
and its Ramifications in September 1998. 



Acknowledgements 



My PhD candidature at The University of Queensland over the period 1995-1998 was 
jointly supervised by Mark D Gould and Jon R Links. In 1992 these two published 
the original description of the Links-Gould invariant, however its evaluation and 
hence the experimental deduction of its properties, remained infeasible. 

In 1997, Louis H Kauffman of The University of Illinois at Chicago visited our group, 
and introduced us to state models for link invariants. This theory, developed by him 
(amongst others) as an elegant method of evaluation for the Jones polynomial, was 
also appropriate for the Links-Gould invariant. 

The confluence of these two streams led to the opportunity for me to explore the 
properties of the Links-Gould invariant. 



Abstract 



The "Links-Gould invariant" is a two-variable Laurent polynomial invariant 
of oriented (1, 1) tangles, which is derived from the representation of the braid 
generator a associated with the one-parameter family of representations with 
highest weights (0,0 | a) of the quantum superalgebra i75[5rZ(2|l)]. We use an 
abstract tensor state model to evaluate the invariant, as per the construction 
of the bracket polynomial state model used by Louis Kauffman to derive the 
Jones polynomial. This model facilitates both computation and theoretical 
exploration. 

Evaluation of the invariant for any particular link involves first drawing a 
presentation of the link as a (1, 1) tangle (open diagram) in a quasi-Morse 
function form. To this diagram is associated an abstract tensor expression 
that is a contraction over all the free indices of a product of rank 2 and 
4 tensors describing the structure of the link. This expression represents a 
summation over all states in a state space, and the primary (rank 4) tensor 
involved is a representation of u. 

In principle, the quantum R matrix associated with any finite-dimensional 
irreducible representation of any quantum (supcr)algcbra will suffice to yield 
an appropriate representation of a. After obtaining a particular instance of 
such a tensor, the model also requires a set of four rank 2 tensors (i.e. genuine 
matrices), and these are chosen using considerations of graphical consistency 
and the properties of quantum (super)algebras. 

The family of four dimensional highest weight [/g[g'Z(2|l)] representations la- 
belled (0, I a) has a two- variable quantum R matrix (unique up to orthogonal 
transformations). Choosing this R matrix to yield the representation of a en- 
sures that our polynomial invariant will also have two variables, and hence 
might be expected to have a greater discriminatory power than an invari- 
ant with only one. We construct this R matrix from first principles, invoking 
Mathematica to assist with the algebra. 

We have cvahiatcd the invariant for several critical link examples and numer- 
ous other links of special forms. Performing the contraction over the abstract 
tensor involves quite lengthy symbolic computation, and is infeasible for a 
human to do. Again, the assistance of Mathematica has been invoked. 

Prom the experimental data, we observe that the Links-Gould invariant is 
distinct from the two-variable HOMPLY polynomial in that it detects the 
chirality of some links where the HOMPLY fails. As we have no method of 
automatically computing the invariant for large numbers of links, we cannot 
make empirical observations about how often our invariant will be able to 
distinguish chirality. Notably, it does not distinguish inverses, which is not 
surprising as we are able to demonstrate that no invariant of this type should 
be able to distinguish between inverses. It also does not distinguish between 
mutants. 
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1 Overview 



The announcement in 1985 of the discovery of the Jones polynomial p4| prompted 
vigorous research which has since yielded many other link invariants. In particular, 
the quantum algebras as defined by Drinfel'd |jl6| and Jimbo being examples 



of quasi-triangular Hopf algebras, provide a systematic means of solving the Yang- 
Baxter equation, and in turn may be employed to construct representations of the 
braid group. From each of these representations a prescription exists to compute 
invariants of oriented knots and links [|^,^,^. The simplest example of this process 
is that the Jones polynomial is recoverable using the simplest quantum algebra 
Uq[sl{2)] in its minimal (2- dimensional) representation. 

From such a large class of available invariants it is natural to ask if generalisa- 
tions exist, with the view to gaining a classification. One possibility is to look to 
multiparametric extensions in order to see which invariants occur as special cases. 
A notable example is the HOMFLYQ invariant which includes both the Jones 
and Alexander-Conway invariants [|],|10| as particular cases as well as the invari- 
ants arising from minimal representations of Uq[sl{n)] Another is the Kauffman 
polynomial which includes the Jones invariant as well as those obtained from the 
quantum algebras Uq[o{n)] and Uq[sp{2n)] in the g-deformations of the defining rep- 



resentations 49 . 



The work of Turaev and Reshetikhin shows that the algebraic properties of 
quantum algebras are such that an extension of this method to produce invariants of 
oriented tangles is permissible. (A tangle diagram is a link diagram with free ends.) 
An associated invariant takes the form of a tensor operator acting on a product 
of vector spaces. Zhang [B3 has extended this formalism to the case of quantum 
super algebras, which are Z2 graded generalisations of quantum algebras. 

Since quantum superalgebras give rise to nontrivial one-parameter families of ir- 
reducible representations, it is possible to utilise them for the construction of two- 
variable invariants. This was first shown by Links and Gould [^] for the simplest case 
using the family of four dimensional representations of f/g[5'/(2|l)]. It was also made 
known that a one-variable reduction of this invariant coincides with a one-variable 
reduction of the Kauffman polynomial by the use of the Birman-Wenzl-Murakami 
algebra. Extensions to more general representations of quantum superalgebras are 
discussed in [EOl. 



^ The HOMFLY polynomial is named by the conjmiction of the initials of six of its discoverers 
, omitting those ( "P" and "T" ) of two independent discoverers |Q . Przytycki, the omitted "P" , 
has furthered the entymological spirit with the suggestion "FLYPMOTH" p256] , which includes 
all the discoverers and has a muted reference to the "flyping" operation of the Tait, Kirkwood 
and Little - the original compilers of knot tables. A similar possibility is "HOMFLYPT" . Bar- 
Natan (Prasolov and Sossinsky p36] cite Bar-Natan Q , who cites "L Rudulph" ) goes further, 
adding a "U" for good measure, to account for any unknown discoverers, yielding the unpalatable 
"LYMPHTOFU"! 



1 



The original paper describing the Links-Gould invariant [3C] does not contain eval- 
uations of it, for want of an efficient method of computation. As such, none of 
its properties were known. Here, we describe a 'state space' method of evaluation 
and investigations into the properties of the invariant. We provide many examples 
and a complete description of the state model for the invariant in abstract ten- 
sor form. This description of the invariant directly facilitates the construction of a 
computer program in Mathematica for evaluation of the invariant. We empha- 
sise that evaluation of the invariant is computationally expensive, and as such it 
would be infeasible to do manually for any but the simplest knots. The code is fairly 
straightforward, and to some extent can be understood directly, without prior ex- 
perience with Mathematica. The readability is somewhat reduced in places where 
computational efficiency has been improved by more sophisticated constructs. 



The method of evaluating the invariant involves a prior construction of the quan- 
tum R matrix associated with a family of four dimensional representations, and 



this is done in explicit detail in § [4.2| . We have used Mathematica to check the 
details. (As an aside, we further illustrate the power of Mathematica by also ex- 
plicitly constructing the associated trigonometric R matrix.) Having obtained this 
R matrix, the construction of the invariant follows from properties of ribbon Hopf 
(super)algebras and their representations. Here we consider the invariants of (1, 1) 
tangles for the following reason: for invariants derived from representations of quan- 
tum superalgebras with zero g-superdimension, the corresponding invariant is also 
zero. If the representation is irreducible, the quantum superalgebra symmetry of the 
procedure ensures that the invariant of (1, 1) tangles takes the form of some scalar 
multiple of the identity matrix. (See |^ for a discussion of this symmetry.) We take 
this scalar to be the invariant. 



We prove that the Links-Gould invariant is not able to distinguish a knot from its 
inverse (Proposition^, nor can it distinguish between mutant links (Proposition]^). 
However it does distinguish many links from their reflections (see Propositions [l| and 
and it is distinguished from the Kauffman and HOMFLY (and hence the Jones) 
polynomials by this behaviour (see § [7.4| and § [7.5| for specific examples) . Whether it 
is a complete invariant for chirality is an open question. 



2 Link Invariants from Combinatorial Topology 



2.1 Basic Terminology of Knot Theory 

A good introduction to the basic concepts of knot theory is found in a book by 
Adams [Q, and aspects of the theory that are particularly pertinent to this work may 
be found in books by Kauffman p6|,^. Between them, these sources describe the 
relationship between physical and mathematical knots, and the concepts of regular 
and ambient isotopy. Here, we assume knowledge of these basics. 

To be certain, multicomponent knots are called links, and almost all of what we 
shall have to say under the heading of 'knot theory' actually applies to links. We 
are in general interested in the classification and properties of oriented links, except 
that we only investigate the invertibility properties of single-component links (i.e. 
true knots). It goes without apology that we are only considering tame (i.e. non- 
pathological) links. 

We shall write "=" to denote the ambient isotopy of link diagrams, meaning that 
they are equivalent under all of the Reidemeister moves ^ ; and we shall use 
the Alexander-Briggs notation and ordering for prime knots (i.e. knots not able to 
be decomposed into a 'product' of simpler knots, which are called composite) . 

Substantial lists of data (i.e. diagrams, values of invariants) for all prime knots 
of up to 10 crossings are found in a book by Kawauchi ||3^, §F]; this information is 



complemented by data for braid presentations in a paper by Jones p5|, pp381-388]. A 
smaller amount of similar data for some (multicomponent) 'prime links' is contained 
in a book by Adams ^ and a paper by Doll and Hoste |jl5[ . 



2.2 The Link Classification Problem 



Consider an arbitrary presentation (i.e. a link diagram) of an arbitrary link. Under 
ambient isotopy, the diagram will in general correspond to a composition of standard 
link presentations, from which we can name the link. Even in the simplest case, 
where the link is prime, how do we discover which standard presentation? This is a 
difficult problem when tackled with only the machinery of the Reidemeister moves, 
and is only made worse when non-prime links are permitted. Thus, whilst we do 
have a classification system for prime links (and a way of decomposing composite 
into prime links), we do not have an efficient method for determining, from a given 
link diagram, to which of the standard presentations our link corresponds. 

Link (polynomial) invariants attempt to address this deficiency. To each link is 
assigned an algebraic expression (typically a Laurent polynomial), and the method 
of evaluation of this expression ensures that it is invariant under the Reidemeister 
moves, i.e. that it is independent of the presentation. 



However, invariants may in general assign the same polynomial to distinct links, 
which reduces their utility. A major goal of knot theory is then to develop a complete 
invariant that assigns distinct polynomials to distinct links. Currently, this goal is 
unsatisfied. The best-known and perhaps the most powerful polynomial invariants 
currently known are the two-variable HOMFLY and Kauffman polynomials, but 
they are far from complete. 

For a particular link invariant to be complete, it is minimally necessary that it must 
distinguish links under the transformations of refiection, inversion (for oriented links) 
and mutation. 



2.3 Reflection, Inversion and Mutation of Links 

Denote by K* the reflection (or mirror image) of a link K. A link is chiral if it is 
distinct from its refiection: K* ^ K; i.e. there are actually two distinct links with 
the same name. Alternatively, a link is amphichiral if it is ambient isotopic to its 
refiection, i.e. K* = K. Note that this definition doesn't require an orientation. To 
illustrate, 3i (the Trefoil Knot, see Figure |r^) is chiral: (3i)* 7^ 3i. Perusal of the 
table of all prime knots of up to ten crossings shows that most of these knots are 
chiral, and it is probably safe to claim that most links are chiral. 

The HOMFLY and Kauffman polynomials can distinguish many (but not all) chiral 
knots from their reflections. The flrst chiral knot that neither the HOMFLY nor the 
Kauffman polynomial can distinguish is 942, i.e. (942)* 7^ 942, but the polynomials 
are equal. Similarly, IO48 is chiral, but the HOMFLY polynomial fails to detect this. 



although the Kauffman does detect it p9|, p218]. (Note that the diagram for IO48 is 



wrongly labeled IO79 in that volume. IO79 is in fact not chiral.) 

Denote by the inverse of an oriented (true) knot K, obtained by reversing its 
orientation. Whilst this is a simple concept for a true knot, there are of course many 
possibilities for the reversal of only some components of oriented, multi-component 
links; we shall not go into these here. Commonly, K = K~^, and we say that K is 
invertible. For example, the Trefoil Knot is invertible (3i)~^ = 3i. Less commonly, 
K 7^ K~^, and we say that K is noninvertible. The flrst example of a noninvert- 
ible prime knot is 817. To date, no known invariant of Vassiliev type can detect 
invertibility. Our invariant is of this type. 

Both the reflection and the inverse are automorphisms of order two, i.e. {K*)* = K 
and {K~^) ^ = K. The notions combined yield: (i^*)"^ = [K'^)* . For example, for 
3i, we have two equivalence classes: 3i = (3i)~"'^ and (3i)* = ((3i)~^) = ((3i)*) ^ . 

A third transformation of a link (diagram) is by mutation, in which a component 
is isolated and replaced with a 180° rotation of itself. In general, this process yields 
a new link. The HOMFLY and Kauffman polynomials always fail to distinguish 
between mutant links. 



A 



2.4 The Links— Gould Invariant 



The paper by Links and Gould |36| describes the construction of hnk invariants from 
a representation of the braid generator a, corresponding to any particular finite- 
dimensional irreducible representation of any particular quantum (super) algebra. 
That paper demonstrates the existence of the Links-Gould invariant, but does not 
supply an efficient method for its evaluation. Here, we set up a state space model 
for this purpose. This model uses an explicit representation of cr, and also requires 
'cap' and 'cup' matrices fi^ and 15^. Below, we develop the state space model for 
the general case of oriented links, before returning to the specific choices for the 



Links-Gould invariant. It is emphasised that we do not reiterate the material of |36 
here, as it is not germane to the state space model. 



2.5 State Space Models 

In order to construct a state space model for a given link, we must first construct an 
abstract tensor expression corresponding to that link, which consists of a contraction 
(i.e. a sum) over the indices of a product of tensors that reflects the link's structure. 
This may be regarded summation of states. 

In order that the contraction yields a topological invariant of ambient isotopy, the 
tensors included within the abstract tensor must satisfy a set of minimal graphical 
consistency considerations; these constraints are well-documented in a paper by 
Hennings p2|, pp59-60] . That paper describes the situation for unoriented links, but 



the material carries over naturally to the oriented case (see a book by Kauffman ||2^, 
pp 235-237]). It is well-known that the use of a representation of the braid generator 
will suffice for the most important of these tensors. 

Further introductions to state space models include the expositions by Kauffman, 
e.g. [^, pp395-399] and ||28|, pp204-209], which evaluate invariants by graphical 



considerations, and ||HT|, §§4.1-4.2, ppl63-172], which evaluates the Jones polynomial 
using the minimal Uq[sl{2)] representation, using only abstract tensors. 



2.6 Morse and Quasi-Morse Functions 



For a link diagram, a Morse function applied to that diagram is a drawing of it in 
terms of vertically-aligned crossings (i.e. with all arrows pointing upwards), vertical 
lines and horizontal arcs p2[. This assignment is not unique. We shall here define a 
^quasi-Morse' function as one that does essentially the same thing, but allows non- 
vertically-aligned crossings and non-straight vertical lines. Our link diagrams (see 
Figures |l^ to ^Oj) are drawn according to such functions; state space invariants are 
quickly seen to be independent of the distinction between Morse and quasi-Morse 
functions. 



2.7 Construction of Abstract Tensors 



To define the abstract tensor corresponding to a particular link, we begin with a 
presentation of the link's diagram according to a (quasi-) Morse function, and assign 
to the arms of each crossing in the diagram a set of labels a, b, c, d, etc. 

By a ^positive oriented^ or 'right-handed' crossing within the link diagram, we intend 
a crossing such that if the right hand thumb points in the direction of one of the 
arrows, the (curled!) fingers of the right hand will point in the direction of the other 
arrow. The opposite situation is naturally called a ^negative oriented^ or 'left-handed' 
crossing. If the arms of a positive oriented crossing with upward-pointing arrows 
have labels a, b, c, d, we shall associate a rank 4 tensor a with that crossing, where 
the positioning of the indices in the tensor corresponds to those on the diagram. 
We shall do the same using cr~^ for negative oriented crossings. Similarly, to each 
horizontally-oriented arc, we shall assign a rank 2 tensor (i.e. a true matrix) called 
a 'cap' or 'cup' as appropriate. This information is summarised in Figure |1]. 




a b a b 




Figure 1: The tensors a and a and cap and cup matrices Q and U , and associated 
oriented crossings and horizontal arcs. 

The abstract tensor for our link is then formed as a contraction over the indices 
of the product of the tensors representing the components of the diagram. This 
amounts to a summation over all possible states, where each state amounts to a 
particular assignment of values to the dummy indices in the terms of the abstract 
tensor expression prior to contraction. To be a little more precise, the construction 
of the bracket polynomial (and hence the Jones polynomial) involves a contraction 
over all the indices of the tensor, whilst our invariant will require a contraction over 
all but two of the indices. This corresponds to drawing the link as a (1, 1) tangle, 
with free ends at the top (index y) and bottom (index x), to which we will assign 
equal, constant values (see § |2.8|) . 

Each crossing in the link diagram will then have a copy of a rank 4 tensor a or 
associated to it, in a manner describing the entanglement of the link. The choice of 
a suitable a is not unique, but a representation of the braid generator associated 
with a quantum R matrix (i.e. a solution of the quantum Yang-Baxter equation) 
will suffice. Evaluation of the abstract tensor will require that a explicitly, and we 
will also require the set of matrices and U^. The latter will be determined from 
other properties of the representation used to represent a. 



2.8 Forming the Invariant from the Abstract Tensor 



Having formed the abstract tensor {Tk)^. for any particular link (where the indices 
X and y correspond to the lower and upper loose ends of the tangle), we form an 
invariant by setting x and y to be the same: 



LG 



K 



K 



(no sum on i), 



(1) 



for any allowable index i; we shall typically choose i = 1. This situation is in contrast 
to the situation for the bracket polynomial, where the invariant is formed as a trace 
(i.e. by summation on i in (|ID). In our case, that trace would have to be replaced 
with a sn^^ertrace as we are dealing with a quantum s«j»eralgebra, which would lead 
to an uninteresting invariant that would always be zero. Figure ^ demonstrates how 
the closure of a (1,1) tangle (open diagram) to form a proper link relates to the 
abstract tensors. 




K I 



iTK)l ■ (n-), , ■ (0+) 



Figure 2: Closure of a (1, 1) tangle to form a closed link diagram. 



In we define the Links-Gould invariant by ([l|), where we have made specific 
choices for the tensors composing Tk- Such an invariant is typically a Laurent poly- 
nomial in one or two variables; indeed our invariant is, although we supply no proof. 



Recall the choice of a representation of the braid generator for a. Deduction of 
appropriate cap and cup matrices follows from graphical consistency considerations 
which ensure that our invariant is automatically an invariant of ambient isotopy, i.e. 
it is writhe-normalised cf. the case of the bracket polynomial pl02]. Figure |^ 
depicts removal of a single twist from a diagram, showing that this requires that we 
must choose a, and 13^ such that: 



,b c 



5^. 



(2) 



We will return to this result in 55.2 




Figure 3: Equivalence of a single loop with an upward line. 
2.9 Auxiliary Tensors 

For any particular link, we may write down its corresponding abstract tensor expres- 
sion using only the tensors cr and cr^^ (as well as Qr^ and O^). However, subsequent 
evaluation of the abstract tensor typically involves the repetition of computations 
reflecting repeated patterns of crossings in the link. To improve computational ef- 
ficiency, we define auxiliary tensors. In the following, we use the notation X to 
represent either a or a^^, corresponding to a diagram component of a crossing with 
upward pointing arrows. The Einstein summation convention is used throughout. 

• The first auxiliary tensors represent crossings that have been 'twisted' relative 
to X. The left, right, and upside-down-twisted versions of X are called X/, Xr 
and Xd respectively (see Figures ^ and ^. 
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• The second set of auxiliary tensors represent copies of the same crossing X 
atop one another (see Figure We use the shorthand 

{x^)i: = {x):).{x^-%i, N^2. 




n 




a c 




Figure 6: The auxiliary tensors X^, expressed in terms of X and X^~^; X is either 
a or (7^^ and N ^ 2.li all arrows are reversed, the definition also holds for X being 
aa or a^^. 



The third set of auxihary tensors correspond to the situation where a crossing 
X (upright, so we write = X for consistency) sits to the left or right of its 
own 'upside-downness' X^ (see Figure |^). 
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Figure 7: The auxihary tensors XdXu and XuXa] X is either a or a ^. 



The fourth set of auxihary tensors describe the situation where a crossing Xi 
is placed atop or below a crossing Xr (see Figure |^). 

[^l^rJb d - y^Oe f ■ y^r)b d 

The diagram for XiX^ is of course a right rotation of the diagram for XdX^- 

The fifth set of auxiliary tensors describe the situation where a chain of 
crossings is formed by the recursive placing of X alongside Xd, with X as the 
leftmost and rightmost crossings; where X is either a or a^^. The recursive 
definition of such chains X^^^ is provided in Figure 
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Figure 8: The auxiliary tensors X/X^ and X^Xi] X is either a or a ^ . The diagram 
for X^Xi is obtained by interchanging r and I. 
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Figure 9: Recursive definition of the chains X^^^ used in the presentation of the pret- 
zel knots; X is either a or a~^. The minimum is the case N = 1, which corresponds 
to X, i.e. XI,, 4 X. 
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3 A Collection of Especially Interesting Links 



3.1 Several Interesting Ordinary Links 



We have selected several links as test examples for evaluation of our invariant. The 
first few examples illustrate the principles involved in the construction of the abstract 
tensors, and the need for the use of the auxiliary tensors. The latter three are 
important cases that we must investigate to determine how our invariant compares 
with other invariants. In Table |l|, we list the links to be studied, in order of their 
Alexander-Briggs notation The data for proper (i.e. single component) knots 
have been taken from §F.2] and for two-component links from (itself citing ||46 



and 15 



K 


Amphichiral? 


Invertible? 


Oi (Unknot) 


Yes 


Yes 


2l (Hopf Link) 


No 


Yes 


3i (Trefoil Knot) 


No 


Yes 


4i (Figure Eight Knot) 


Yes 


Yes 


5i (Cinquefoil Knot) 


No 


Yes 


52 


No 


Yes 


(Whitehead Link) 


No H, pp49-50] 


Yes 


6i 


No 


Yes 


62 


No 


Yes 


63 


Yes 


Yes 


7i (Septfoil Knot) 


No 


Yes 


72 


No 


Yes 


817 


Yes 


No 


942 


No 


Yes 


IO48 


No 


Yes 



Table 1: Data for the links to be investigated, see Figures to 
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Below, we describe each of our example links in terms of (1, l)-tangle (open diagram) 
forms, where the indices x and y correspond to the lower and upper loose ends of the 
tangle. From these diagrams, their abstract tensors are immediately obvious. 

Some of the diagrams have been constructed using a braid presentation, as taken 



from |2^, ppl09-110] and |2^, pp381-386]. For the Unknot (Oi), a braid presentation 



is the trivial e G -Bi as the tangle is a simple vertical line, and we have (TqJ^ = 



2\ (Hopf Link) and 3i (Trefoil Knot): Diagrams, constructed from the braid 
presentations cxi^ and cxi^ G -B2, are found in Figure |TD|. The abstract tensors 
follow immediately: 



X 



b c 



b c 





Figure 10: Tangle forms of 2\ (the Hopf Link) and 3i (the Trefoil Knot). 



4i (Figure Eight Knot): From the diagram found in Figure |TI], we have: 



d e 



5i and 7i (Cinquefoil Knot and Septfoil Knot): Diagrams, constructed from 
the braid presentations (j\ G B2 and cr[ G i?2, are found in Figure |T2|. We have: 



{"11 



y a 



7\y a 

b 



(IS-) 



b c 



b c 



1 A 




Figure 11: Tangle form of 4i (the Figure Eight Knot). 



A diagram is found in Figure |T3|. We have: 
5^ (Whitehead Link): This hnk is named after the topologist J H C Whitehead 



p200] . A diagram is found in Figure [T^. To deduce the abstract tensor, we 
initially define a temporary tensor to reduce computation: 



With this, we have: 



d b 



6i: A diagram is found in Figure |T5|. We define an auxiliary tensor: 
With this, we have: 



1 c 




Figure 12: Tangle forms of 5i and 7i (the Cinquefoil Knot and the Septfoil Knot). 




Figure 13: Tangle form of 52- 




Figure 14: Tangle form of 5f (the Whitehead Link). 



62: A diagram is found in Figure |T5|. We define an auxihary tensor: 
With this, we have: 



63: A diagram, drawn from the braid presentation cr 2^1 cr2 G -B3, is found in 
Figure |16|. We define some auxihary tensors: 



With these, we have: 



, k c 



1 ^7 




1 o 




1 n 




l-i- A diagram is found in Figure |T^. We have: 



817: A braid presentation is [o^^o^^ cr2'^r^^2 £ -B3, from which we draw the dia- 
gram found in Figure 18. Again, we define some temporary (rank 6) tensors 
to reduce computation: 



With these, we have: 

To reduce computation, we may define even more auxihary tensors: 
where: 



k I \ Jm n 



on 



[EA) 



[EB) 




Figure 18: Tangle form of 817. 



01 



: A diagram is found in Figure 19. Again, we define a temporary tensor to reduce 
computation: 



With this, we have: 




Figure 19: Tangle form of 942. 



IO48: A braid presentation is cr^ ^cr|cr^ '^a2 G -B3, from which we may draw the 



diagram found in Figure Again, we improve efficiency by the use of some 
auxihary tensors: 



(TA) 
(TB) 



o- y 9 

b d h 



b d h 



)b f-{'' )dh 

-3\b d / \e h 



Then: 



1048; 



(TB) 



b d h 



i k 



00 



3.2 A Class of Noninvertible Pretzel Knots 



The first demonstration of tlie existence of noninvertible knots p25] was by 
Trotter in 1964 who showed that the well-known class of pretzel knots labelled 
{p,q,r) in the notation of Conway |T^, were noninvertible when p,q and r are 
distinct, odd, and greater than 1 (their ordering is irrelevant). Trotter is of the 
opinion that the knots are chiral (for some evidence, see § |7.5| ). This class of pretzels 
provides an easily-programmable set of examples for us to evaluate our invariant to 
see if it detects their noninvertibility. 



The structure of the knots (p, g, r) in this family is depicted in Figure pTl where the 
auxiliary tensor is used. 




p Ni y 

udu/b e 



I Udu^ y"ndu)c / 






I I 

I ^Udu' y"udu)d X 
I I 

X 




Figure 21: Tangle form of the (noninvertible) pretzel knots {p,q,r) 



The abstract tensor associated with the pretzel is: 



or; 



\a y _ 



udu'b 



/ q \b e 
y'^uduJc f 



c f 
uduJd X 



O A 



3.3 The Kinoshita— Terasaka Pair of Mutant Knots 



The 'Kinoshita-Terasaka Pair' is the best known example of a pair of mutant knots 
that are known to be distinct ( [|I], pl06,174] and |^). To be precise, the first of 
the pair is usually known as the "Kinoshita-Terasaka Knot", and the second has 
been called the "Conway Knot" , as Conway used it in coining the term "mutant" . 
In the original source by Kinoshita and Terasaka @, pl51], the knot involved is the 
one labelled k{2, 2) (reproduced in |3^, p53]). They had constructed this knot as an 
example of a nontrivial 11 crossing knot with Alexander polynomial equal to 1. The 
source used to draw our example is from [|l], pl74]; note that these diagrams have 
12 crossings, so they are not minimal. 

It is known that neither the HOMFLY nor the Kauffman polynomial can distinguish 
any pair of mutants [|l|, pl74], but stronger statements can be made when these in- 
variants are regarded as being examples of Vassiliev invariants. A Vassiliev invariant 
is defined p229] to be of type d if it is zero on any link diagram of ci + 1 nodes, 
and of degree d if it is of type d but not of type d — 1. In 1994, Chmutov, Duzhin 
and Lando 0] proved that all Vassiliev invariants of type at most 8 will agree on 
any pair of mutants. The Links-Gould invariant is certainly a Vassiliev invariant, 
although we do not know its type, but it seems unlikely that it is of type greater 
than 8. We shall investigate whether the Links-Gould invariant distinguishes these 
mutants. 



We illustrate KT, the first of Kinoshita-Terasaka pair, in Figure ^ where the 
components KTA and KTB are defined below, in Figures ^ and From KT^ we 
may build the mutant KT' by replacing the component KTA with KTA' (depicted 
in Figure ^), which is formed by reflection of KTA about a horizontal line. 




Figure 22: KT^ the first of the Kinoshita-Terasaka pair of mutant knots, where 
the subdiagrams KTA and KTC are found in Figures and ^ respectively. The 
mutant KT' of KT is obtained by exchanging KTA with KTA' (see pSf). 



The tensors associated with KT and KT' are: 

{Tkt)1 ^ {KTA)l\.{KTCti-{n-)^^.{U^)-.{n^)^^.{l5^Y 

{Tkt')1 ^ iKTA'):l{KTC)ti-{n-),,-{u^r'-{n^)^^-{u^Y 

where: 

{KTBtl 4 (a)f^(a^t^„■K-V-):;■(^^+),,■(0-)^- 



{KTB) 




Figure 23: The component KTC of the Kinoshita-Terasaka pair. 




Figure 25: The component KTA' of KT\ the second of the Kinoshita-Terasaka pair. 
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4 The Uq[gl{2\l)] Representation (0, | a) 



We consider the family of four dimensional representations of highest weight (0, | a) 
of the quantum superalgebra ?7q[5fZ(2|l)] (for nonzero complex q 'generic', i.e. not a 
root of unity). These representations depend on a free complex parameter a, and 
are irreducible for a ^ 0, —1, which we shall assume throughout. Below, we use the 
notation of which itself is based on the nonquantised case gl (2|1) found in 0. 



Material in this section has been obtained by specialising that in §0. 



The material in this section is presented in great detail, with many routine deriva- 
tions presented in full. This is to ensure correctness, and to illustrate the principles 
used to explicitly construct these quantum superalgebra representations. In §|, we 
discuss the idea of the automation of the construction of quantum superalgebra 
representations. 



4.1 Some Standard Notation 



The following operations are actually defined for any quantum superalgebra Uqlg] 
(where g is a Lie superalgebra), in particular, for Uq[gl{m\n)] (see 



4.1.1 The Graded Commutator 



We shall employ the graded commutator [■, ■] : Uq[gl{2\l)] x Uq[gl{2\l)] Uq[gl{2\l)], 
defined for homogeneous x,y E Uq[gl{2\\)\ by: 

[x,y] = xy- {-)^''^^y^yx, (4) 
and extended by linearity. 



4.1.2 The q bracket 



We introduce the q bracket, a shorthand notation defined for x G [/5[5'/(2|l)] by: 



A q"" -q 



The q bracket should not be confused with the graded commutator bracket! Observe: 
lim = X. 



on 



4.2 The Quantum Superalgebra Uq[gl{2\l)] 
4.2.1 Z2 Grading on f/J^/(2|l)] Indices 

Firstly, we define a Z2 grading [■] : {1, 2, 3} Z2 on tlie tliree gl (2|1) indices: 

[1] = [2] = (even) 
[3] = 1 (odd). 

Tliroughout, we shall use dummy indices a, 6 = 1, 2, 3 where meaningful. 



4.2.2 Uq[gl{2\l)] Generators 



We shall follow the notation used by Zhang ppl237-1238] to define Uq[gl{m\n)], 
except that we shall substitute^ the notation q^~^ for Ka- Specialising to the 
case Ug[gl{2\l)], our superalgebra has 7 simple generators: 

qE\^qEh^qEh (Cartan) 
E'^2,E'^i (raising) 
i?^3, E^2- (lowering) 

Note that we are not writing the Cartan generators as -E"a- The fact that there are 
2 simple raising generators indicates that ?7q[5fZ(2|l)] has rank 2. 

In addition to the above 7 simple generators we have 2 others defined in terms of 
them H, pl971, (3)] and ||, pl238, (2)]: 



(a) ^ E\E^;-q-^E\E\ (raising) 

(6) E'\ = E\E\-qE^iE^2 (lowering) 



(6) 



On all of the 9 generators we define a natural Z2 grading in terms of the grading 
on the indices: 



[E\] = [a] + [b] (mod 2), 

where the definition of the exponential shows that we have [q^""] = [E°'a] = 0. The 
product of homogeneous x,y E Uq[gl{2\l)] has degree: 

[xy] = [x] + [y] (mod 2) . 



^After this substitution, our notation for the Cartan generators approaches that of another of 
Zhang's papers ppl970-1971], which hsts E°-a as the Cartan generators. Strictly speaking, this 
causes problems with the meaning of 9^ " as a formal power series. 



on 



4.2.3 U^[gl{2\l)] Relations 



For the simple Uq[gl{2\iy\ generators, we have the following f/q[(7/(2|l)] relations: 



The Cartan generators all commute: 

^ q±E\^±E^.^ q^'^q-^"^ = 1. (7) 



The Cartan generators commute with the simple raising and lowering genera- 
tors in the following manner: 

From (I), we have the following useful interchange: 
Replacing q with q~^ in @ yields the equivalent: 

^_(_)W^.,^,^^^ = (10) 

The squares of all odd generators are zero (this will also apply to E^i and £"^3, 
although we will never actually use this in practice): 

{E\Y={E%Y = 0. (11) 

The non-Cartan generators satisfy the following commutation relations (this 
is the really interesting part!): 



[E^2,E'^i\ = [E^i—E'^2 

[E\ES] = [E\ + E%] 

[E\,ES] = 

[E\,E\] = 0. 



q 



(12) 



There are two additional relations, known as the Serre relations. Their inclu- 
sion ensures that the algebra is reduced enough to be simple. We include them 
for completeness, but they will not concern us further. 



{E\yE\- {q + q-')E\E\E\ + E\{E\Y = 
{E\YE'2-{q + q~')E\E'2E\ + E',{E\Y = 0. 

These relations may be more succinctly expressed using the two nonsimple 
generators, viz: E\E\ - q^E^sE\ = and - q-^E^iE^i = 0. 

More generally, expressions involving the two nonsimple generators E^i and E^s 
may be manipulated by prior expansion using @. 



Q1 



4.2.4 Ug[gl{2\l)] Simple Generator Interchange Rules 

Application of the definition of the graded commutator (|[) to the LHS of (|l^) yields: 

Substitution of all possible choices for a and h in (|T3|) yields the following explicit 
commutation expansions: 



[E\,E\] = E\E\-E\E\ 

[£^^3, -E^2] = E'^^E'^2 + E^2E'^3 

[E^2-,E^2] = E^2E'^2 — E^2E^2 

!rn2 rp2 ] _ rp2 rp2 rp2 rp2 



(14) 



Combining the information of the commutators of ([T^) with their expansions in ( JT^ ) 
yields the following interchanges ('commutations'): 



[a 



E\E\ = [E\-E\]^ + E\E\ 



(b) E\E\ = [E\ + E\\-E\E\ 

(c) E\E\ = E\E\ 

(d) E'^sE'^i = E'^iE'^3. 



Lastly, expansion of @ and (|T0|) to all possible cases yields: 



(15) 



(a) 




= q^E^q'^'^ 


q~ 


-^'^E\ 


= q-^E\q-^'^ 


(b) 


q^'^E\ 


= E\q^'^ 


q' 




= E\q-^'^ 


(c) 


q''"^E\ 


= q-^E^q"""' 


q' 


'^'-E\ 


= q'E\q~''"' 


id) 


q^"'E^2 


= E\q^"^ 


q~ 


-E'-E\ 


= ESq-^"' 


(e) 


q^"'E\ 


= q-^E\q^^^ 


q~ 


-^'-E\ 


= q^E\q~^^^ 


(/) 


q'''^E\ 


= q'E\q^'- 


q~ 




= q-'E\q-^'^ 


ig) 


q^'^E\ 


= q'E\q^'^ 


q' 




= q-'E\q-'''^ 


(h) 


q'^'^ES 


= q-^ESq"""' 


q~ 




= q^ESq-'''^ 


(0 


q'^'^E^ 


= E\q^'^ 


q~ 




= E\q-^'^ 


(j) 


q^'^E\ 


= q~^E\q^'^ 


q~ 


-EhE\ 


= q^E\q-^'^ 


(k) 


q^''^E\ 


= E\q^'^ 


q~ 


-EhE\ 


= E\q-^'^ 


(0 


q^'-^ES 


= q^ESq^'^ 


q~ 


-E'^ES 


= q~^ESq-^'^ 



} (16) 



The definition of the nonsimple generators (^, together with the Cartan commu- 
tations d^), the nilpotency in (|Tl|) and the exchanges in (pTSj) and (0), provide all 
the necessary information for manipulation of vectors encountered in the explicit 
construction of representations for Ug[gl{2\l)]. 



4.3 The Uq[gl{2\l)] Module V 



We shall be concerned with the ?7g[(7/(2|l)] module V = V(o,o|o)5 which is irreducible 
and typical (recall that we insist a ^ 0, —1). We may obtain a concrete realisation of 
it using the Kac induced module construction (from which we immediately discover 
that ^ is 4 dimensional). This section follows [|19|; with the caveat that there is a 
small error in that work, and that we shall choose a slightly different definition for 
our basis. Where differences occur, we shall mention these. 



4.3.1 Dual Bases and for V 

Let denote a basis for the four dimensional Uq[gl{2\l)] module V. In accor- 

dance with the Kac induced module construction, we shall define these basis vectors 
as created by the action of all possible nonzero combinations of all (not just the 
simple) the odd lowering generators on defined to be the highest weight state: 

(a) |2) 4 P2E\-\1) ] 

(b) |3) ^ f3sE\.\l) \ (17) 

(c) |4) 4 (3,ESE\.\1), J 

where the = 2,3,4 are suitable normalisation constants, to be deduced be- 
low (implicitly, /?i = 1). We shall use the symbol "■" to represent algebra-module 
multiplication. 

The above definition (P!7| ) differs from that (implicitly) used in in that the 
definitions of |2) and |3) are interchanged. The choice here is more natural in that 
the vectors are ordered in terms of decreasing weight. To be certain: 

|1) has weight (0, | a ) 

|2) ( 0, -1 |a + l) 

|3) (-1, 0|a + l) 

|4) (_i,_i|« + 2). 



Consistent with the Z2 grading on Uq[gl{2\l)], following []TU|, pl58], we grade the 
basis states by: 

[|1)] = [|4)] 4 0, [|2)] = [|3)] 4 1. 

By "consistent", we mean that the action of a homogeneous [/g[5'/(2|l)] generator x 
on a (homogeneous) basis vector yields a homogeneous vector of degree [x] + [ji)]: 



QQ 



Thus, declaring |1) as even means that |2) and |3), which are created from the 
action of odd lowering generators on |1), must necessarily be odd. Lastly, |4) may 
be regarded as being defined in terms of the action of an odd lowering generator on 
the odd 1 2), hence is necessarily even. 

We define a basis dual to where [(^l] = [\i)]. Written in component 

form, these dual vectors are represented by the transpose complex conjugates of the 



original basis: {i\ = \i) = . Then we have: 

(^I|J) = (^|J) = S.r 



4.3.2 Action of the Generators on the Basis 

The module V has highest weight vector of weight (0, | a). Using the shorthand 
A = (0, I a), we have Ai = A2 = and A3 = a. 

• As 1 1) is a weight vector, the action of the Cartan generators on |1) is:[| 

E\-\l) 4 A,|l), (19) 

i.e.: 

E\-\1) = 0, E\-\1) = 0, E%-\l) = a\l). (20) 

Expansion of the exponential, and the use of (|19D yields: 
g^-|l) = g^>|l). 

We take a moment to justify this, using the observation that, recursively for 
all integers j ^ 2 we have {E\y -11) = {E\Y'^E\ = \i{E\)^'^ ■ hence 
(E^)^^|1) = (A.)^|1). 

^ gin(,)i.% ■ ii) = f; . |i) = Y^^^^{E\y ■ |i) 

Thus, (^) is equivalent to: 
9^^^-|l) = g^^^-|l) = g°|l) = |l), g^^«-|l)=g11). (21) 

■^The equations of ( pO| ) are of course the relations that would define a gl (2|1) representation as 
having weight (0, | a). As we are regarding as the generators of Uq[gl{2\l)], perhaps it 

would be more honest to label this the (f , 1 1 q^°') representation, as per 



OA 



Temporarily writing E^q = 1 and Aq = 1, for any complex scalars 7^: 



1 



■|l) = g^-^-|l), 2 = 0,. ..,3, 



which may be extended to: 

. |1) = (22) 



Using (^) , we have: 



= [ELo7iAi]Jl). 
In our case, where A = (0, | a), this simplifies to: 

[Elo7.^^],-|l) = [7o + 73«],|l). (23) 

• As 1 1) is a highest weight vector, it is annihilated by the action of the ?7g[(7/(2|l)] 
raising generators: 

E\-\l) = E\-\l) = E's-\l) = 0. (24) 

• The even subalgebra of gl{2\l) is gl{2\l)Q = gl{2) (Bgl{l), and the Kac induced 
module construction builds a gl{2\l) module by acting on a gl {2\1)-^ module 
with all the odd lowering generators. This naturally extends to the f/q[(7/(2|l)] 
case. 

As |1) spans a one dimensional Uq[gl{2\l)]Q = Uq [gl (2)] © Uq [gl (1)] module, 
it is annihilated by the action of the even Uq [gl (2)] lowering generator E'^i. 

E\-\1) = 0. (25) 

The definitions of the duals imply that (|l^), (0) and (pSj) have natural analogues: 

(1| ■ E\ = A,(l|, (1| ■E\ = {1\- E\ = {l\.E\ = 0, (1| ■ E\ = 0. 
We shall use these implicitly, without further comment. 
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4.3.3 The Normalisation Constants (3i 



We demand that the basis is normalised, viz: 

= 1, (26) 

and this defines for i = 2,3,4. Let us assume for now that a and q are real 
and positive. (Our results for the R matrix will hold later by analytic continuation.) 
Initially, we declare that |1) is already normalised, i.e. (1|1) = 1, equivalently /3i = 1. 
Before proceeding to deduce the normalisation constants, we deduce some auxiliary 
results, which will also be of use when we come to discover the representation itself. 

[E\ + E\]^E\E\ -11) 

{q - q-T' {q'^'-Q^'^E^ - q-^'^q-^'^E\) E^ ■ |1) 
[q - q-y (q^'-E\q^''-^ - q-^'^ E\q-^''-^) E^ • |1) 

E\{q - q-y' [qVS'''' ' q-^-^^'q-^'^) E\ ■ |1) 
E\{q-q-')-' {^q'q^'-^E^q^'^ -q-'q-^'^E^q-^'^) • |1) 

E\EMq-q-y' (g^g^^^g^^^ - g^^g-^^^g"^^^) • |1) 
= E\E', [1 + E\ + E\]^ ■ |1) @ [a + 1\E\eS ■ (27) 



[E\-E\]ES.\l) 

''^ (g - g-^)^' [q''\~'''^E\ - q-'\''^E\) ■ |1) 

(g - g-^)"' (g^g^^^i^;^2g-^'^ - g-^g-^^^i5^'2g^'^) ■ |i) 

(g - q'^y^ [q^ESq^^'^q-^"^ - g-^S^^g-^'^g^'^) • |1) 
= E\[l + E\ - E\]^- |1) 9 [1]^ES ■ |1) = E^ ■ (28) 

E\ . |1) {E\E\ - qE\ES) ■ |1) ® -qE\E\ ■ (29) 



E\E\.\l) -qE\E\E\.\l) 

-q[[E\-E\]^ + E\E\)E\.\l) 

-g {[E\ - E\]ES + E\ESE\) . |1) 

® -q[E\ - E\]E\ . |1) @ -qE\ ■ (30) 



E\E\.\l) E',{ESE\-qE\E\).\l) 

® -qE\E\ES.\l)^^ -qE\E^,E\.\l) 
-qE\[[E\ + E\]^-ESE\)-\l) 

® -qE\ [E\ + E%1. |1) # -q[a]^E\ ■ |1) # 0. (31) 



^ 3^ 2- |1) 



(156) 



P3h 



([^^2 + E^]^- E\E\) . |1) ® [E\ + E^]^ . |1) 
HJl)- (32) 



{e\[E\ - E\]E\ + E^E23i5'2i?S) ■ |1) 
E\[E\-E\] E\-\l) 



E\E\.\l)^[a]\l). (33) 



We now compute the normalisation coefficients: 

(2|2) = (2| ■ |2) ■ E\E\ . |1) ® [a]M. 



Demanding that (2|2) = 1 then amounts to selecting (32 as real and equal to [a]^ 
(up to a phase factor, which we shall set as 1). 



(3|3) = {?>\-\?>)^^ M;{1\-E\E\.\1) 

®^ q'PMA ■ E\E\E\ES . |1) ® qMMl 



Demanding that (3|3) = 1 then amounts to selecting /^s as real and equal to 
(again, up to a phase factor, which we shall set as 1). 



on 



(4|4) 

= (4| ■ |4) ®^ /34/?4*(l| ■ E\E^,E\E\ ■ |1) 

/34/34*(l| ■ ^^^3 {[E\ + E% - ^3,^23) E\ ■ |1) 
® g2/34/54*(l| ■ E^E\ {[E\ + E% - E\E\) E\ES ■ |1) 

q'm(^\ ■ {e\E\[E\ + E\\E\ES - E\E\E\E\E\E\) • |1) 

g'/JA* ■ i?^3i?^2[i?'2 + ES]^E\E\ . |i) - ■ Ei^i^^ • |i)) 

® g2/54/34*(l| ■ i5'3i5'2[i?'2 + ES\E\E\ ■ |1) 
® g2[a + l\m{^ ■ E\E\E\ES ■ |1) 
® g^NJa + l]/4/54*- 

Thus, up to a phase factor, we have: 13/^ = + l]^"^''^- In conclusion, we 

have the normahsation constants: 

(&) /^s = q-M',"'" 
that is, the definitions of ( P^D may be replaced by: 

W |2) ^ N;^/^i?^2-|i) 

(6) |3) ^ ^ (35) 




QO 



4.4 The Representation tt = ttx 



In terms of these dual bases and {(^|}^^]^, we define a representation tt = ttx 

on Uq[gl{2\l)]. For the generators, we have: 



(a) 


r^{E\) 
vr {E\) 


= -|3)(3|-|4)(4| 


ib) 


= -|2)(2|-|4)(4| 


(c) 


vr (^33) 


= + + l 


id) 


7^{E\) 

rr{E\) 


= -|2)(3| 


(e) 


= -|3)(2| 


(/) 


vr iE\) 


= Hf|l)(2| + [a 


(9) 


7r{ES) 


= Hf|2)(l| + [a 



|2)(2| + (a + l)|3)(3| + (a + 2) |4)(4| 



l]f |3)(4| 
l]f |4)(3|. 



>(36) 



This comes from |T^, pl57, (3)], although that paper does not describe any of the 
conventions concerning the definition of the basis vectors, or the phase factors in their 
normalisation constants. The interchange of the definitions of |2) and |3) between 
that paper and this work generates the same changes in (|36|) . (Our paper |jl3|, p3] 



uses the conventions of [DJfl rather than the above. 



Observe that we have used the expression E°-a here, rather than the more technically 
correct g^"" (or even g*^")'"'^""). The representation is of course extended to the rest 
of [/g[5'/(2|l)] by linearity. In particular, for future reference, we have: 



(a) vr(g^^) 

(b) 7r(g^^^) 

(c) niq^''-^) 

(d) 7t{E\) 

(e) n{E\) 



|l)(l| + |2)(2|+g-i|3)(3|+g-M4)(4| 
|l)(l|+g-M2)(2| + |3)(3|+g-M4)(4| 
g° |1)(1| + g"+i |2)(2| + |3)(3| + 
-[« + l];/^|2)(4|+g-M<^|l)(3| 
l]f |4)(2|+gHf |3)(1|. 



|4)(4| 



} (37) 



\a 



The expressions presented in 
where X e Ug[gl{2\l)]: 



may be constructed using the following definition. 



vr(X) 4 X-vr(/)=X-(|l)(l| + |2)(2| + |3)(3| + |4)(4|) 

= {X . |1))(1| + (X ■ |2))(2| + (X ■ |3))(3| + (X ■ |4))(4|. 



(38) 



In § |4.4.1| and § [4.4.2| , we evaluate the action of the ?7q[5'/(2|l)] simple generators on 
the basis vectors, then use (^) to construct the explicit representations presented 
in 



on 



4.4.1 Representation of the Cartan Generators 



Firstly we consider the Cartan generators. InitiaUy, we deduce the form of vr(g^^i^ 



q^'^-\l) ^'^^ (39) 

|2). (40) 

g^'^-|3) Psq'''^E\.\l)9-q'P,q^'^E\E^.\l) 

®^ -P,E\q^'^E% . |1) -P,E\E^q^'^ . |1) 

© -PsE\ES.\l)^9q-^P,E\.\l)^\-^\3). (41) 

g^^^-|4) P,q^'^E-\E\.\l)^^P,E-^q^'^E\.\l) 

® g-/?4/53-^E^,.|3)@'^-Ai^;^,E^.|l)&^-|4). (42) 



Combining the resuhs of (|39|) to (|42|) , we have, using (pS]), in accordance with (^g): 

vr(g^'^) = + |2)(2| + g-M3)(3| + g-M4)(4| , 

equivalently, in accordance with (|5B|a): 

7r(E^) = -|3)(3|-|4)(4|. 

Let us pause to justify this statement. Starting with (P7|a), firstly, we have, as E^i 
is a diagonal operator/matrix: 7r(g^^i) = q^^^ ^\ thus: 

7r(g^^^) = e'-^M^'^) = ^ (-^'i)^ ^ ^ [In (g)]-[7r (E\)]- ^ 

by the definition of the exponential as a formal power series. Next, we have, by 
inspection: 



He': 



)V {E\) if i > 

/ if i = 0. 



A n 



Thus, we have: 



7r(g 



j^ [\niq)nniE\)r ^ ^_ ^ ^^^^^ ^ [In (g)]' 



+ / 



i=l 



1=1 



— vr 



-vr 



i=l 



E 

.i=0 



+ 1 



= -n {E\) [e'-^'^r^ - ij + / = (1 - g-i) n {E\) + I 

= (1 - q-') (- |3)(3| - |4)(4|) + |1)(1| + |2)(2| + |3)(3| + |4)(4| 
= |l)(l| + |2)(2|+g-M3)(3|+g-M4)(4|. 

A similar expansion may be made in the opposite direction. Below, we shall simply 
report results for the other Cartan generators without further ado. 

We repeat the process for the other Cartan generators. For n^q^^^): 





|1) 


(21) 




|2) 








^) 




|3) 


^) 












© 




|4) 











P,q^'^E^ . |1) @) q-%ESq^'^ ■ |1) # q-'p,E^ ■ |1) 

Psq^^'^E-^-ll)^ -q'P,q^'^E\E^.\l) 
-q'(3,E\q^'^ES ■ |1) -q^ (3,E\E\q^'^ ■ |1) 
-q'hE\E\.\l)^ P,E\.\l) |3). 



(43) 
(44) 



(3,q'^'-ESE\ . |1) »^ q-^(3,ESq''''-E\ ■ |1) 



g-i/?4/53-^i?^ ■ |3) g-^/34i^^2£;^ ■ |i) g-^|4). 



Combining the results of ( ^31 ) to (^6]), we have, in accordance with (|3 
vr(g^'^) = + ?-M2)(2| + |3)(3| + g-M4)(4| 



(45) 
(46) 



thus, we have, in accordance with (|36|&) 
7r(i?22) = -|2)(2|-|4)(4|. 



y( 1 



Lastly, for 7r(g^ ^): 







(21) 






•|2) 


^) 


/?2g^'^i?^2 -11)^1 






&) 


g"+i|2). 




■|3) 


&) 








^) 








© 


-q\''(3^E\ES 




■14) 


(0=) 


P,q^'^ESE\ . 1 











(47) 

q'P,E^q^'^.\l)^qVP,ES.\l) 

(48) 

|1) -q'PsE\E^q^'^ . |1) 

)®g"+V^/33i?^-|l)^^"^^|3). (49) 
g^g"+^/?4/?3"'^'2 ■ |3) q^^^P,E\E\ ■ |1) & g"+^|4). (50) 



Combining the results of ( ^7D to (^0]), we have, in accordance with (|3 
7r(g^^«) = |1)(1| + |2)(2| + |3)(3| + |4)(4| , 

thus, in accordance with (p^c): 

TT (E^g) = « |i)(i| + (a + 1) |2)(2| + (a + 1) |3)(3| + (« + 2) |4)(4| 



4.4.2 Representation of the Non-Cartan Generators 

Now, we deduce the forms of the simple raising and lowering generators. Firstly, let 
us deduce the form of tt (£^^2): 



E\- 




(24) 


E^2 ■ 


■|2) 


(§«) 


E\^ 


■|3) 


^) 








E\- 


■|4) 


^) 






(§) 



(51) 

P2E\E\ . |1) & P2E\E\ ■ |1) @ 0. (52) 
P,E\E\ . |1) @ -qP,E\ . |1) -g/33/32-^|2) 
-|2). (53) 
^,E\E\E\.\l)^^ P,E\E\E\.\l) 
-ql3^E\E\-\l)^\. (54) 



Combining the results of ( pID to (^4]), we have, in accordance with (|3 
7r(i?^2) = -|2)(3|. 



As we know that the representation is unitary 



0, 



21 



, we have 



rc{E\)=n{E\y 



where, in matrix notation, X'^ is the conjugate transpose of X. Thus, immediately, 
in accordance with (p6|s): 



{E\)=-\3){2\ 



Secondly, vr (-E^a) is more work: 





|1) 


(24) 


E\ 


|2) 






|3) 






|4) 


&) 






























(06) 



0. (55) 



/3,E^3i5^,.|l)»^«]J/^|l). (56) 



P,E\E\ . |1) ® 0. (57) 

(3^e\eSe\-\i) 

-qP, [[E^ + E^]^- ESE\) E\E\ . |1) 



-g/34 ([a + l],^'i^'2 ■ |1) - E\E\E\eS ■ |1)^ 
-g/54 ([a + l]^E\ES ■ |1) - [a]^E^E\ ■ |1)) 
+ 1]^(3,E\E'2 ■ |1) ® [« + l]/4i?^ ■ |1) 
[« + l]A/53-'|3)^''^[« + l]ri3)- (5^ 



Combining the results of ( p5|) to (^8]), we have, in accordance with (|36|/): 
^(^'3) = Hf |l)(2| + [« + l]y^|3)(4|. 



Again, from the unitarity of the representation, vr (-E'^2) = (i^^s)^, hence, in accor- 
dance with 



^(^'2) = [<'|2)(l| + [a + l];/^|4)(3| 



A Q 



We do not really need to discover the forms of it (E^^) and it (E^i); however we 
mention them here for completeness. In accordance with (^rf, e), we have: 



TC 



TT 



TT 



{E\) TT {E\) - q-\ {E\) TT (^^2) 
|2)(3|-(H;/^|l)(2| + [a + l]f |3)(4| 

+g-i(Hf |l)(2| + [a + l];/^|3)(4|).|2)(3| 
+ |2)(4|+g-iHj/^|l)(3|. 
[E\) TT [E\) - gvr [E\) tt [E\) 
(H;/^|2)(l| + [a + l]f |4)(3|)-|3)(2| 

+g|3)(2|.(Hf |2)(l| + [a + l];/^|4)(3| 
|4)(2|+gHf |3)(1|. 



A A 



4.5 ?7q[5'/(2|l)] as a Hopf Superalgebra 



Ug [(7/(21 1)] may be regarded as a Hopf superalgebra, when equipped with suitable 
coproduct A, counit e and antipode S structures [0, pl971]. Here, we only describe 
the coproduct, as that is all that we shall require below. Details for the counit and 
antipode may be found in ^M. 



4.5.1 Coproduct A 

We define a coproduct (a.k.a. comultiplication) structure, which is a Z2 graded 
algebra homomorphism A : Ug [5-/(211)] Ug [5/(211)] O Ug [5/(2|l)], by: 



a) A{E\) = I®E\ + E\®I, a = 1,2, 3 

b) A{E\) = E\^q--2^''"'-''''^ + qh(^''-''''^ ®E\ 

d) A{E^s) = E%®q-"2(''''+''''^^ + q"2(^''+E'^^ ®E\ 

e) A{E^) = E\®q-"2i^'^+^'-^^ + q^^^^'^+'''^^ ®E\, 



and extended to an algebra homomorphism on all of Ug [5/(2! 1)]. As A is a homo- 
morphism, we have of course that A(l) = 1 ® 1, and that for x,y E Ug [5/(2|l)]: 

A{xy) = A{x)A{y). (60) 

A being graded means that it preserves grading, i.e. [A(x)] = [x] for homogeneous 
X E Ug [5/(211)]. Of course, (^a) implies that: 

A(g^''") = g-^"" (g) g-^"". 



4.5.2 The Twist Map T and the Alternative Coproduct A 

As well as this standard coproduct, there exists another possible coproduct structure: 
A, defined by A = T- A, where T : [/J5/(2|1)](8)[/J5/(2|1)] f/j5/(2|l)]Ot/j5/(2|l)] 
is the twist map, defined for homogeneous elements x,y E Ug[gl{2\l)]: 

T{x®y) = (-)[^][J'] {y®x). 



We shall require this operator below, in § [1.8.1| . More generally, the notion of a twist 



map is applicable to the tensor product of any two vector spaces. 



4.5.3 A{E\) and A{EK^) 



In § 4.6.2 and § 4.6.3 , we shall require the services of A{E^i) and A(i?^3) respectively. 



AiE\] 



(66) 

(g) 
(E|c,e) 



nm,g,h,k) 



A{E\E\ - qE\E^2) 
A{E^)A{E\)-qA{E\)A{E^) 
{e\ ® + q\iE^.+ES) ^ e3 

-q (e\ ® q-^.iEh-Eh) + q^,(Eh-E^,) ^ 

(E^aE^i - q^E\E\) ® q-^^'^+Eh) 
^q\(E\+E^,) ^ I^eSE\ - q^E\E\) 

- (g - q'^) q'2(i^'^~^'^^E^ ® E\q-"2(^"'+^'^^ 
E\ ® q-hi.E\+Eh) + qhi.E\+ES) ^ E\ 

- (g - g-i) g^^^'i-^'^)^^^ ® E\q-"^^^"^+^'^\ 



E^2 



(61) 



(6a) 

i),(i) 



(|li^e,/,i) 



da) 



A{E'2E\ - q~'E\E\) 
A{E\)A{E\)-q-^A{E^,)A{E\) 
{e\ ® + q\iE\-E\) ^ ^1^ 

[e\ q-UE\-E\) ^ q\{E\-E\) ^ ^1 

{E\E',-q-'E',E\)®q-l^'''^^'''^) 
^^liEh+Eh) ^ (^i^^2^ _ q-^E\E\) 

+ {q- q-') q'^(^''+^'-^^E\ ® E^q--^^"^'^-"^"'^ 
E\ ® g-|(^'i+^'3) + ^|(i^^+i^^3) ^ 



(62) 



4.6 The Tensor Product Module V (^V 



From the 4 dimensional representation on the space V , we may construct a 16 
dimensional representation on the tensor product space V ®V . The grading on the 
basis for V induces a natural grading on the basis for V ®V: 

[K) ® m = im + m ■ 

The Z2 graded action of Ug [gl{2\l)] ® ?7g[(yfZ(2|l)] on K ® K is defined for generators 
X, ?/ G f/q[5'/(2|l)] and vectors u,v & V in the following manner: 

{x ® y) ■ {u ® v) = (— )'^^'"^ (x ■ n y ■ t>) , homogeneous y, u, (63) 

and naturally extended by linearity. 

4.6.1 Decomposition of the Tensor Product Module 

The tensor product module V has the following decomposition with respect to 
the coproduct for generic values of a: 

V^V = Vi®V2®V3, (64) 

where: 

Vi has highest weight (0, | 2q; ) and dimension 4 

V2 ( 0,-1 I 2a + 1) 8 

V3 (-1,-1 I 2a + 2) 4. 



This decomposition is found in p2770]; which, strictly speaking, refers to the non- 
quantised case gl (2|1), but which extends naturally to the quantised case Uq[gl{2\l)], 
cf. H, pl91]. 



We shall construct symmetry adapted orthonormal bases {|\E'j)}^^| and 
for the spaces Vi and V3 respectively, in terms of the basis elements of V. 

In the construction process, we shall introduce arbitrary phase factors; these will be 
mentioned in passing. The bases are: 



(a) \^\) = 

(b) m) = + (gf|l)® |2) + g-t|2)®|l)) 

(c) m) = (g° + g-")-^(gt|l)®|3) + g-t|3)®|l)) ^(55) 

(d) =^ (g" + g-")-^[2«+l]-i 

[a + ® |4) + g-"|4) ® |1)) - N|(g5|3) ® |2) - ® |3)) 



«]|(g"+'|4) ® |1) + q-^-'ll) ® |4)) + [a + l]|(g5|3) ® |2) - g-i|2) ® |3)) 
(6) 1^3^ = + ('g^|4)^|2)+g-^|2)® |4)') J>(66) 

(C) 1^3^ = ^g^|4)^|3)+g-^|3)^|4) 

(rf) |vl/l) = |4)®|4). 



Observe that these submodules do not have purely odd or even bases, but that each 
basis vector is homogeneous. 



Our basis is similar to those presented in |T9[ and |jl3[, except that we have made 
some different choices in the definitions of the Of those two papers, the basis 
presented in the former contains a slight error as well. In particular, to get from the 
basis in [jl3[ to the above, we must interchange the definition of |2) with |3) (and 
also |\I'2) with l^'g) and |\E'2) with I^E'D), and replace |4) with — 14). These changes 
naturally lead to some different phases in the respective normalisation constants. 

Bases dual to {|\E'^)}^^i (for k = 1,3), labelled {(^^^1}^^^ respectively, are found 
from the definitions: 



= k = l,3, j = l,...,4, (67) 

= (-)[l^)Hb-)l((,|®(^-|), i,j = l,...,A. (68) 



Observe that the general form of each of the basis vectors l^'^) is: 



where the coefficients O^^m are in general complex scalar functions of q and a. To 
simplify the discussion, we shall limit q and a to be real and positive, which ensures 
that the 6^^ rn are real, and hence equal to their own complex conjugates. In that 



case, from (|67|) and (|68|), the duals of these vectors |\E''j) are given by: 



Explicitly, we have: 



(a) {n = 

(6) = + (gt(l|® (2|+g-t(2|®(l|) 

(c) (^^1 = + (3|+g-t(3|®(l|) 



> (69) 



[a + ® (4| + g-"(4| ® (1|) + (3| ® (2| - g-2 (2| ® (3|) 



(a) (^fl = + g-"-i)-5[2a + 

[a]|(g"+^(4| ® (1| + ® (4|) - [a + l]|(g5(3| ® (2| - g-5(2| ® (3|) 

(6) (^3| _ (^a+i^^-a-i^-U^H±i^4|^^2|+g-^(2|® (4| 



>(70) 



(4|®(4|. 



The multiplication operations between the dual bases for V (E)V are given by: 



(a) (|z)®|j))((A:|®(/|) 

(&) ((^|®(j|)(|fc)®IO) 



^,j,A;,/ = 1,2,3,4. (71) 



So, where do the bases |\Ef*^) come from, and why aren't there basis vectors for V2? 



These things are explained in § [4.6.2| to §4.6.4| . 



A n 



4.6.2 An Orthonormal Basis for Vi 



The highest weight vector of is |1), of weight (0,0 | a), hence the highest weight 
vector of K ® V" is |1) ® |1), of weight (0, | 2a). We shall label by Vi the submodule 
ofV<^V generated by this vector |1) ® |1), i.e. Vi has highest weight (0, | 2a), from 
which we may determine, using Kac's dimension formula, that it has dimension 4. 
We shall construct an orthonormal basis for Vi, using the Kac induced 

module construction. Initially, as already defined in (|65|a): 

The normality of |1) implies the normality of |1) |1): 

® (ii) ■ ® ii)) i-f''''''' ((111) ® (111)) i-ri = 1, 

For Vi, the equivalent of (||) ioi V ig) V is A{E'^i) ■ |^|) = 0. This fact will be 
used implicitly. The other basis vectors for Vi are generated from the action of (the 
coproduct of) the odd lowering generators of Uq[gl{2\l)], in much the same manner 
as the basis for V itself was generated from cf. (PT]): 

(a) m) 4 P^,A{ES).{\1)®\1)) ] 

(b) m) 4 (31A{E\).{\1)^\1)) \ (72) 

(c) ^ (31A{E^E\) . ® \1)) , J 

where the for j = 2, 3, 4, are suitable normalisation constants. As the coproduct 
is a homomorphism (|5D|), we have A {E^2E^i) = A (-E'^2) A {E'^i). Again, the basis 
vectors are ordered in terms of decreasing weight: 

I^J) has weight ( 0, | 2a ) 

1^2) ( 0, -1 I 2a + 1) 

1^3) (-1, 0|2a + l) 

1^1) (-1, -1 I 2a + 2). 

Without further ado, let us discover |\1'2)- Firstly, we have: 

A {E^) . ® \i)) 

& E^ -11)® . |i) + .\i)^eS-\1) 

® E32-|l)®g-t|l) + gt|l)®E=^2-|l) 
[<^/'(g^|l)®|2) + g-t|2)®|l)). 



To deduce /^g, we demand (\l/2|\I/2) = 1- Investigating: 



[gt(i|® (2|+g-f(2|® (1|] ■ [gt|l)® |2) + g-t|2)® |1)] (g° + g^") 



hence we require, up to a phase factor Pi = {q°' + q ^ . Thus, (|72| a) yields, in 
accordance with ('|65|&): 



1^2) = + Q'T' ® |2) + ?"^|2) ® |1)] • 



To build |^'3), we evaluate: 

(El) (E^E\ - q'E\E^) ■ |1) ® g-|(^^i+^-^3) . |i) 



m -q^E\ES . |1) ® g-t |1) - g^+f |1) ® E^E^^ ■ |1) 
- (g - g-i) g^(^^-^^^)E32i5^32 ■ |1) ® ■ |1) 

(r^i^;^ . |2) ® |1) + gf |1) ® E\ . |2)) 
gH|[gt|l)®|3) + g-t|3)®|l)]. (73) 



To deduce we demand (^3 1^3) = 1. Investigating: 

[gt(l|® (3|+g-t(3|® (1|] • [gt |1) ® |3) + g't |3) ® |1)] (g" + g"") 



hence we require, up to a phase factor /3| = (g" + g ^ . Thus, ([72^) yields, in 
accordance with (|65| c): 



1*3) = + 9"")"' [9^11) ® |3) + g~^|3) ® |i)] . 



In building l^^l), we require: 



A{E^)A{E\).{\1)^\1)) 



(Hell 



(MM?) 



gHl (E 



q[a 



q[a 



q[a 



q[a 



(g-^|3) ® |1) ® |3)) 



+gtE32-|l)®g-5(^'2+i?33).|3) 

+gfg5(S'2+i?^3).|l)^^3^.|3) 



gt [a + 1]| 1 1) ® |4) + g-t [a + 1]| |4) ® 1 1) 

-g-fg|(^'2+i?33).|3)^[c,]||2) 

+gt[a]||2)®g-i(^'2+^'^)-|3) 
[a+l]|(gt|l)®|4)+g-f|4)®|l)) 

-N| (g-tgl+f |3) ® |2) - gtg-5-f |2) ® |3) 

[a + l]|(g"|l)®|4) + g-"|4)® |1)) 
-[«]|(gi|3)®|2)-g-5|2)®|3)) 



To deduce Pi, we demand (\l/4|\l/4) = 1. Investigating: 



[a + l]|(g°(l|® (4|+g-"(4|® (1|) 
+ H|(gi(3|®(2|-g-i(2|®(3|) 



[a+l]|(g11)® |4) + g-"|4)®|l)) 
-[a]|(g5|3)®|2)-g-i|2)®|3)) 

[« + 1], + + K + 9"') 

(g - q-y' [g'^"+i - g-^-i + _ 
(g - g-^)"^ (g" + g"") (g2°+^ - g^^^-i) 
[2a + 1], (g- + g-") , 



hence we require, up to a phase factor P} = q ^[a]^ ^ (^q» _)_ g 2 pa + 1]^ 2 . Thus, 
([T^) yields, in accordance with (§^(i): 



[a + 1]| (g"|l) ® |4) + g--|4) ® |1)) - (g^|3) ® |2) - g-^|2) ® |3)^ 



CO 



4.6.3 An Orthonormal Basis for V3 



The construction of a normalised basis for V3 is, in a sense, dual to that of Vi. The 
lowest weight vector of V is |4), of weight (— 1, — 1 | a + 2) (see (0)), hence the 
lowest weight vector of V 1/ is |4) (g) |4), of weight (-2, -2 | 2a + 4). 

We shall regard the submodule ofV^V generated by the action of (the coproducts 
of) the odd raising generators on this vector |4) ® |4) as V3, i.e. V3 has lowest weight 
(—2, —2 I —2a + 4), from which we may determine, again using Kac's dimension 
formula that V3 has dimension 4 (as did Vi). 

We shall construct an orthonormal basis {l^'^)}^^^ for V3, using the Kac induced 
module construction. Initially, as already defined in (|65|a): 

= |4) ® |4). 

As for |1) ® the normality of |4) implies the normality of |4) |4): 
((4| ® (4|) • (|4) ® |4)) = (-)[<^I][I^>1 ((4|4) ® (4|4)) = {-fl = 1. 



The other vectors in V3 are generated from the action of (the coproduct of) the odd 
raising generators of ?7g[5'/(2|l)] on this lowest weight state = |4) |4): 



(a) 


m) 


^ PIA{E\).{\A)0\A)) 






4 /5|A(Ei3) ■ (|4) ® |4)) 


(c) 


m) 


^ f3!A{E\E\).{\A)®\A)), 



(74) 



where the for j = 3, 2, 1 are suitable normalisation constants. Again, the vectors 
are ordered in terms of decreasing weight: 

has weight 



-1, -1 I 2a + 2) 
-1, -2 I 2a + 3) 
-2, -1 I 2a + 3) 
-2, -2 I 2a + 4). 



Recall that |4) is a lowest weight vector of Vi. This means that it is annihilated by 
all lowering generators. Without further ado, let us discover l^'l). Firstly, we have: 



A{E\) -(14)® |4)) 



f) 
(§/) 



E% . |4) ® g-|(^'2+^'3) . |4) + q^iEh+Eh) . |4) ^ E\ ■ |4) 



[a+l]f |3)®g- 



\iE\+Eh) 



[a + l]y' [g'^K) ® |3) + g-'^|3) ® |4) 



To deduce we demand (\l/3|\l/3) = 1. Investigating: 



^"^(41 ® (3| +g-'^(3| ® (4| ■ g^|4) ® |3) + g-'^|3) ® |4) 
(@>) ^ ^ 



hence we require, up to a phase factor /3| = [a + 1]^^'^ {q°'^^ + q ° ^) Thus, 
(fr^g) yields, in accordance with (|6^c): 



To build I^E'f), we evaluate: 



A -(14)® |4)) 

(E^E^s - q-'E\E\) . |4) ® g-|(^^i+^^3) . |4) 
+qliE\+E^) . |4) ^ (^i^^2^ _ g-iE^gi^;!^) . |4) 
+ {q- q-') qhiE'-2+Eh)E\ . |4) ® E^q-'^^''"'-^"^^ ■ |4) 



-[a + l]]'^ g"*"|4) ® |2) +g~"*"|2) ® |4) 



(75) 



To deduce we demand (^E'il^l) = 1- Investigating: 



Q + 1 



g— (4| ® (2| (2| ® (4| ■ g— 14) ® |2) + g"— 12) ® |4) 



hence we require, up to a phase factor /3| = [a + l]^^^(g"'*"^ + g " ^) Thus, ([74|6) 
yields, in accordance with 



1^3) = + g^|4)® |2) + g-^|2)® |4) 
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In building l^l/f), we require: 



A{E\E\) -(14)® |4)) 



(ii/§,c) 



A{E\)A{E's) -(14)® |4)) 



1/2 



E\ 



^i("+i)|4) ® |2) + g-5(°+i)|2) ® |4) 



-[a + 1] 



1] 



1/2 



1/2 



gi(°+i)E23 . |4) ® g-5(^^'2+s«3) . |2) 

+g|{°+l)g|(^'2+^'3). |4)^^2^. |2) 

+g-5("+i)E23 . |2) ® q-¥.E'^+Eh) . |4) 

[(AT [g"+'|4)® + |4)] 

+ + U|3)®|2)-g-i|2)®|3)' 



To deduce we demand (\I'f |\I''^) = 1. Investigating: 



[g"+'(4|® + (4|] 
-[a + l]f [gi(3|®(2|-g-i(2|®(3| 

,1/2 



Hr [g"+'|4)®|i) + g-"-^|i)® |4)] 

+ [a + 1];/' [g5 13) ® |2) - g-^ |2) ® |3) 

N, (g^"^^ + g-^"-^) + [a + 1], (g^ + g-^) 



(g° - g"") (g'°+' + g~'"~') 

+ (g°+i - g-"-i) (gi + g-^^ 



(g - g-')"' [g'""'' - g"'""' + g" - g"1 = (g"+' + g"""') [2« + 1],. 



Thus, up to a phase factor (5\ = + g " -'^) ""^^^^ 
yields, in accordance with (|66| a): 



1/2 



I*?) = (g"^' + g 



■[2« + i]; 



1 

2 . 



HI (g^-^'K) ® |1) + g-"-^|l) ® |4)) - [a + 1]| (g-^|2) ® |3) - g^|3) ® |2) 



4.6.4 No Basis for 1^2? 



Although it is not needed below, it would be more complete to find such a basis. V2 
is necessarily 8 dimensional, and its highest weight is (0, — 1 1 2a + 1), determined as 
the highest remaining weight after the weights used in the basis for Vi have been 
eliminated. In order to build a basis for V2, we would have to define a highest weight 
vector \ml) as a linear combination of weight vectors of this weight, and then build 
the rest of the basis by the action of (the coproduct of) the lowering generators on 
\ml). In fact, we may write down l^'^) immediately. The only possible combinations 
of weights that sum to (0, —1 1 2q; + 1) are those of the vectors |1) ® |2) and |2) (g) 
This combination already appears once, in {mD; examination of that vector (|65|&) 
shows that the only other possible orthogonal, normalised combination is (as usual, 
up to a phase factor): 

1^2^ = (g" + g--)-5 (g-f |1) ® |2) - gt|2) ® |1)) . 

In this case, the underlying gl (2|1)q module is 2, not 1 dimensional, i.e. the action 
of the even lowering generator E'^i (viz A{E'^i)) on l^'^^ is not to annihilate it, cf. 
(p5|). Thus, the construction of V2 would involve requiring the next lowest weight 
vector of the top level to be: 

\ml) 4 PiA{E\).\ml). 

This vector is actually of a lower weight (viz (—1, | 2a + 1)) than a vector in the 
next level, hence its numbering is peculiar. This would complete a basis for the top 
level of V2. 

The next level would have a basis {1^1/2), 1^4), I ^5); l^e)} consisting of the 4 vectors 
formed by the action of E^2 and E'^i on \ml) and where our vectors are 

naturally ordered in terms of decreasing weight. So far, we have: 



\ml) of weight 


( 0, 


-1 


2a + 1) 


m) = f3lA{ES) . \ml) 
\ml) = f3iA{E\) . \ml) 
\ml) = /3lAiE\)-\ml) 


( 0, 


-2 


2a + 2) 


(-1, 





2a + 1) 


(-1, 


-1 


2a + 2) 


\n) = f3lA{E^).\ml) 


(-1, 


-1 


2a + 2) 


\ml) = (3lA{E\).\ml) 


(-2, 





2a + 2) 



Observe that the natural definition of V2 in this manner yields a basis that contains 
a pair of vectors of the same weight, viz l^^^) and 1^1). This pair will not naturally 
be orthogonal, so we may have to use the Gram-Schmidt process to orthonormalise 
them (although we might be lucky!). 

The construction of a basis for V2 would then continue for one more level. 



4.7 Projectors onto the Subspaces C F 1^ 



From the sets of basis vectors and their duals we may construct 

projectors Pk V ®V ^ Vk by: 



Pk = J2\^j)^^% ^ = 1'3' (76) 

where the multiphcations are defined by (^). Of course, we could also define P2, but 
we haven't defined a basis for V2 . . . The gentle reader is reminded that projectors 
satisfy: PiPj = SijPj. 

Now let / be the identity operator on the tensor product space V tSiV, defined as: 



As Pi + P2 + P3 = I, we thus do not need to explicitly construct P2 (or even a basis 
for V2); we simply set: 



P2 ^ I-Pi- F3. (77) 

An alternate notation is to let e*j = be an elementary rank 2 tensor, which 

ensures that we must write the elementary rank 4 tensors as: 



e'^,i = e',®e\ = \i){j\®\k){l\. (78) 

We then have: 

4 4 

/=$^e^®e^,= 5^e^^,,. (79) 



We may substitute the definitions of the dual bases for V ®V (|6^, ^ and |70|) 
into the definition of the projectors Pi and P3 ([76|) , paying heed to the definitions of 
the tensor multiplications (|7T|a) to yield the following explicit expressions (checked 
with Mathematica!): 



Pi 



+(g" + g" 



fe^\2 + g-"e2i2i) 



+ (g°ei3^3 + g-^e^igi) 



+ (g" + g-")-'[2« + ir- 
r + ((g2V^4 + g-'Vl4l) + (e^^4 + e^l)) 1 



M|[« + i]|- 



2 6 23 

-a+i 32 
2 6 41 

g,a+2g32^^ 



+ N, {{q'e^\2 + g-^e2323) - {e^\, + 6^33^)) 



> . 



(80) 



+ q 

q 



[2a + 1]-' ■ 



N, ((g'"+'e«4i + g-'"-'eiS4) + ie'\, + e'\i)) 



/ g"+2e^^23-9~ 



-i 14 



_^q,a+2g32^^ — q- 



'2 623^4 



+g " 2 6^^14 — g"+ 2 6^3^,^ 
\ +q 23 — q 2 6 32 y 



+ [« + 1]^ ((gie3232 + g-ie2323) - (632^3 + 6^332)) 



-(g"+i + g-"-i) 
(g"+ie42^2 + q- 



-a-l„24 \ I /„24 „42 ^ 

e 24 J + (,e 42 — e 24 J 

,34 



I /„o+l„43 I „-a-l„34 \ _i_ /„34 „43 N 

+ [q e 43 + g e 34j + (^e 43 — e 34; 



(81) 
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4.8 The Quantum R Matrix and the Braid Generator 



4.8.1 The Quantum Yang— Baxter Equation 



For any classical Lie superalgebra g, the associated quantum superalgebra Ug [g] 
is an example of a quasitriangular Hopf superalgebra, which necessarily admits a 
universal R matrix R EUq [0]®f/<j [g] satisfying (among other relations) the quantum 
Yang-Baxter equation: 

-R12-R13-R23 = -R23-R13-R12, (82) 

where the subscripts refer to the (standard) embedding of R acting on the triple 
tensor product space Uq[Q\ ® Uq[Q\ ® Uq[Q\. 

From any representation of f/g[0], we may obtain a tensor solution of (|8^) by replac- 
ing the superalgebra elements with their matrix representatives. Observe that the 
R 'matrix' is in fact a rank 4 tensor. Similarly to (|7lD , multiplication of the tensor 
product of matrices a, 6, c, d is governed by: 

(a ®h) {c® d) = (— )[^1W (ac ® hd) , homogeneous 6, c. 

Recall that we were careful to only define our basis for V ®V for q and a real and 
positive. As it is known that the R matrix is unique,^ analytic continuation makes 
our R matrix valid for any complex q and a. 



4.8.2 The Graded Permutation Operator P 



We introduce the graded permutation operator P : V ® V V ® V ^ defined for 
graded basis vectors v'^.v^ G V by: 

P{v''®v') = {- f^^'\v' ®v''), (83) 

and extended by linearity. (We use the shorthand [v^] = [k].) Specialising to the 
case of the f/g[(7/(2|l)] representation tt = tta, in terms of the decomposition (|^) of 
® V^, we findQ that as g — 1, P is the identity operator I on Vi and V3, and — / 
on V2: 



P 



lim [Pi 



Po + P. 



We will require P in §4.8.3 



The R matrix is unique up to several considerations, see p236] or [ p3[ . 
^This is taken from [Q, but with a shght change of notation and a convenient choice of nor- 
malisation. 



4.8.3 Obtaining the R Matrix R and the Braid Generator a 



The quantum R matrix R mentioned in §[4.8. 1| is a specialisation of a trigonometric 



R matrix R{u), which corresponds to the associated affine quantum superalgebra 
f/q[g'^^)] not the basic Uglg], where the parameter u in R(u) is introduced in the 
process of the extension of Ug[Q] to f/g[0*-^'']. As any representation of the affine 
extension is a representation of the underlying quantum superalgebra, R{u) may 
also be regarded as acting on V ® V, indeed R{u) satisfies the spectral parameter- 
dependent Yang-Baxter equation [|19| (cf. (p2D): 



R{u)^^R{u + v)^^R{v)^^ = R{v)^^R{u + v)^^R{u)^^. (85) 

Setting R{u) = PR{u), in it was shown that in our case (i.e. g = gl{2\l)), R{u) 
has the form: 



R{u) oc ^ + ^2 + , P3. 

^ ' ^ qU+2a qU g2a+2 



In [19|, the constant of proportionality was chosen as 1. Here we choose —1; which 



will ensure (see below) that ( p^ is satisfied. Thus: 

^(-) = -rr^^l - - ,._ ,2..2 ^3. (86) 

In passing from f/g[0*^-'^^] to Ug[g], we take the limit as u — > 00, and our R matrix 
does likewise, viz: 



R = lim R{u) = q-^'Pi - P2 + g'^+^Pg. (87) 

u^oo 

Now, we define PB|,P7|: 



a ^ R. 



Taking lim^^oo in the spectral parameter-dependent Yang-Baxter equation (^): 
(a ® /) (/ ® a) (a ® /) = (/ ® a) (a ® /) (/ ® a) . 



It follows that 0" is a representation of the braid generator, which is a sufficient 
requirement that a may be chosen as the (positive) crossing matrix when computing 
our link invariant. 



Using (173) and ®, ( 



p7D becomes: 



= g"" {q" + g"") Pi + + g'""') ^3 - J] e'^^,-. (88) 

ij 

From the explicit expressions for the projectors (l80|j8l|) , we may now construct 
an exphcit expression for a. To simphfy its presentation, we use the substitution 
p — q°', which systematically replaces a, so that we may regard a as a function of 
q and p rather than q and a. We also set: 



Y 



[p ^ - q^ +p^q^ - if^ = [q{p-p ^) {pq - p ^q ^)] 



lMl/2 



1 _ 1 
q2 — q 2 



^9) 



With the help of Mathematica, we find: 



a = p-V^i + (p-^ - l)e% - 

+ ip-' - l)e'\i + iq' - l)e'\2 - e'%s 

+Y\'\, + {p\' - l)e'\, + {p\' - l)e'\,+pW'^ 

+P~\e'\2 - e'\^) + p-\e'\, - e'h{) 

He'\A + e'\{)-q{e'\, + e'\2) 
+Y{e'\,-e'\,) + qY{e^\,-e'\,) 

+pq{e^\2 - e'\,) + pq{e'\z - e''u). (90) 



Inspection of the basis vectors for V ® V shows that we may expect a maxi- 
mum of 36 nonzero entries in a. In fact, a has only 26 nonzero entries; ten entries 
(fortuitously) cancel to zero. 
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4.8.4 Aside: Evaluation of R{u) 



Having constructed projectors, it is a simple matter to apply (^) to them to yield 
the trigonometric R matrix R{u), of significance in the study of exactly solvable 
models 



IS]. Again, we let Mathematica do the work: 



R{u) = -[a-|]Ja + |];V\i-(e2222 + e3^3) 
-[a + l + fya + l-f];W 

-[«] Ja + + f];^[a + 1 - f];^(g"e^i4i + q-^'e'S,) 

+n| [« + 1]| [f] j« + f + « - f];' X 

. (62^41 - e^i^a) + (e^24i - 6^132) 

+ [f] Jl - f] J« + + 1 - firie^i + e''i4) 

+ (g-g-^)"^[a + f];^[l + a-f];^x 

-[a + l]Ja + 1 - f];^ [gt(e42^2 + 6^343) + q-He'\, - e'\d] 
+ [f]'[« + + l-f]r(e''32 + e^3) 

-[f ] Ja + 1 - f ];'(e^\2 - £^^24 + e'\, - e'h,). (91) 



This tensor has 36 nonzero entries: the fortuitous cancellations that caused a to have 
only 26 nonzero entries do not occur here. Observe that lim^j^oo R{u) = R = a. From 
this expression, one may obtain the Hamiltonian H^{U) defined in ||19|| . 
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4.8.5 A Skein Relation Satisfied by a 



From (88), it is straightforward to deduce tliat a satisfies the polynomial identity: 



^- V3 ^ (^-1 _ ^-2.-1 _ ^2a+l) ^2 ^ _ ^-2.-1 _ ^2«+l) ^ + = Q, (92) 

known as a skein relation. To justify (^), begin by observing that from (P?]), the 
eigenvalues of a are —1 and g^""*"^, hence: 

(a - {a + I){a- q'^^'l) = 0, 

which may be expanded to yield: 

+ (1 - - g2a+2^) ^2 ^ (^2 _ ^-2. _ ^2a+2^^ a + I = 0, 

which may be immediately seen to be equivalent to (P^). 

We shall evaluate our link invariant using the state model based on the components 
of a. The invariant would also be able to be evaluated in some cases using (p^), 
although not in all cases since it is of third order. (The invariant may also be directly 
evaluated for a class of links using quantum super algebra representation theoretic 
results 0.) 



4.8.6 Construction of R Itself? 



We have defined a = PR, shortcutting the actual finding of R. For good measure, 
we might obtain R itself from a in the following manner. Observe that where |\E'^(g)) 



is a basis for V as per ( |65| , |66| ) and P is as per (p3|), we have: 



thus where: 



k=i j 



for the eigenvalues Pi = q ^",P2 = — 1,P3 = q'^°'~^^, we have: 



R 



k j k j 

E/^^Ei^?(^"'))<^?(^)i- 
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5 Explicit Construction oi a, a , Q and O 



5.1 The Origin of a 



The explicit form of the representation of the braid generator a on the graded 
module V ^ V was presented in (pO]). The grading on V means that this tensor 
actually satisfies a graded version of the quantum Yang-Baxter equation (|82D , which 
looks similar, except that it includes strings of parity factors. Our abstract tensor 
formalism for the construction of link invariants also requires tensor multiplications 
that do not include such parity factors. Thus, for our purposes, we shall require a 
version of the above a, as it appears after the grading on V is removed. We shall refer 
to that tensor here as a, however it is still labelled "a" as far as the link invariant 
material is concerned . . . 



5.1.1 Components of Rank 4 Graded Tensors 



Let A represent an arbitrary graded rank 4 tensor acting on V ^V, then in terms 



of scalar coefficients A* ^ : 



ijkl ijkl 

Note that we are representing the scalar coefficients as A* ^ rather than the more 
standard f, this is for consistency with the way that the abstract tensors for 
the link diagrams are presented. This does not apply to the elementary tensors. 

Let us write {v''} as the (standard) basis for V, where of course the v'' satisfy: 

e>^ = Sk,v\ (94) 
hence the action of A on V ®V is: 

ikmn 

© V- .i k f^lj](lk] + [n]) (i i I 



ikmn 
ikmn 



ik 



5.1.2 Aside: Components of P 

Specialising (|95D to the case where the operator A is P, we have: 



ik 



ik 



Equating components: P] f = so that 



This expression simplifies; the only nonzero entries are: 



From this we may reconstruct P: 

ijfci ijkl ij 

- 13 , „14 I ^21 22 23 , „24 

— e 11 — e 21 — e 31 + e 41 + e 12 — e 22 — e 32 + e 42 

, ^31 ^32 33 , ^34 , 41 42 43 , 44 

+e 13 — e 23 ~ e 33 + e 43 + e 14 — e 24 — e 34 + e 44. 

A check on the consistency of our results is that from (|88D, we expect lim^ 
Observing that lim^^i y = 0, we see that indeed this is the case. 



5.1.3 Components After Removal of the Grading on V 



In (^), the parity factor is constructed from the degrees of vectors; in the ungraded 
case, there would be no such factor, indeed we would have: 



Note that we are expressly not dropping the grading on A, although we shall refer 
to the new operator as A. This motivates us to set: 

(97) 

Removal of the grading on V has the net effect of changing the sign of one in four 
of the entries in the respective tensors. 

5.1.4 Explicit Components of a 

Applying (|97D to the explicit expression for a (i.e. (pO|)) yields the ungraded edition 
of cr: 

a = p-V^i + (p-^ - l)e^l - e^, 

+p-\e\2 + e'\) + p~\e'\s + 
+ (e^^4 + e^l)-g(e^3 + e^2) 
-r(e^3 + e%)+gF(e^l + e^2) 
+pg(e% + e^^24) + M(e^3 + e^=^34), 

which is the tensor that we shall use to represent the positive crossings in the abstract 
tensor formalism. 
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5.1.5 Explicit Components of a ^ 

Similarly to (^), we may find cr^^ from the knowledge of the eigenvalues of a. Where 
the eigenvalues of a are —1 and g^""*"^, the eigenvalues of are g^", —1 and 
g,-2o-2 respectively, thus: 

To get to the version of after the grading on the indices has been removed, 
we again invoke (pT]). This yields: 

-e'\2 + {q-' - l)e'\, + {p^\-' - l)e'\, 

+ (p- V - l)e^4 + V V\4 
+P(e^\2 + e^i)+p(e^^3 + e^i) 

-q-'Y{e'S2 + e'\,) + q~'Y{e'S, + e'\,) 
+p-\-\e'\, + e'\2) + P'\-\e'K^, + e'\,), 

which is the tensor that we shall use to represent the negative crossings in the 
abstract tensor formalism. Perhaps not surprisingly, also has 26 nonzero com- 
ponents. 

5.1.6 Writing a and as Matrices 

It may be easier on the eyes to view a and in the more 'graphical' form of 
matrices rather than via lists of components. Our rank 4 tensors have 4 indices, 
hence we cannot present them elegantly on the page in matrix form as we can do 
with rank 2 tensors. We can however, present a rank 4 tensor on the two dimensional 
medium of a page as an ordinary matrix using certain standard conventions. (This 
really amounts to an isomorphism.) 

To begin with, we shall represent rank 2 tensors as matrices, that is, the elementary 
rank 2 tensor e*j will be represented by the elementary matrix ejj-. Recall that our 
underlying carrier space is 4 dimensional. This means that all indices will run 
from 1, . . . , 4, i.e. we intend ejj to be a 4 x 4 matrix. 

We adopt the (usual) convention that the elementary rank 4 tensor e^^ ji = e*^ ® e'^i 
is constructed by insertion of a copy of the elementary rank 2 tensor e''i at each 
location of e*j (i.e. each element of e*j is multiplied by the whole of e'^i). This means 
that e^'^ji is represented by the elementary (16 x 16) matrix e4(j_i)+fc^4(j_i)+/. 
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Our convention then tells us that A* ^ appears as the (4(i — l) + k, 4(j — 1) + 
of ip{A) (the matrix representing A), explicitly: 



^(^)4(i-l)+fe,40-l)+Z - ^) I- 

Using (p8|), we evaluate the expressions for the matrices ip(a) and ipi'^^^)'- 



p-2 . . . 


p-1 . . . 


p-i . . . 


1 . . . 


. p-1 . . 


. . . 

. -1 . . 


-q . . 


-Y . . . 

pq 


. . p-K 


. -q . 


. . . 

. q^-1 . . 
. -1 . 


qY . . . 

pq 




. -Y . 
■ ■ PI 


. qY . . 

■ pq 


y2 . . . 

. p^q^-1 . 

. p^q^-1 . 

■ P'^Q^_ 



. . . g-2r2 


p . . 
. . q-^Y 


p . . 

.-q-^Y . 




. p 

. . . g-2y 


.-1 . 
. . <7-2-l 

. . . p-2g-2-l 






. . p 

. . . -q-^Y 




-1 

. . . p-2g-2-l 


p^lg-l 




• • • p"^g"^ 







5.2 Caps and Cups and O 



We will use: 



0+ = h 



-2{a + l) 



-2(c+l) 



„2{q+1) 



2(q + 1) 



These choices for and are sufficient although not unique. To justify them, we 
begin with the comment that the graphical consistency considerations of Figure ^ 
(invariance under the first Reidemeister move) require: 



This means that fl^ = (O^) ^, so we need only find 




X 




±\a y 



5t 



Figure 26: Invariance under the ffist Reidemeister move. 



5.2.1 Choosing U+ and n+ 



may be chosen from graphical consistency considerations involving the quarter- 
turn rotations in Figure ^. The choice: 

i.e. O"*" = /4, ensures that the definition (|), i.e. 



(where X is either o" or o" simplifies to the elegant form: (X^)^ ^ = {Xy^ ^, and 
indeed we also have {Xd)'^ ^ = (X)^ ^. Clearly, we then have = I4, viz 

, = SI (99) 

5.2.2 Choosing ^ and fi- 



For the choice of 15 , we begin by imposing the graphical consistency requirement 
of Figure |^ (invariance under ambient isotopy), i.e. (|): 



abc 



The choice of in ( p9|) simplifies this to: 

E (^)' t ■ (^")' = (100) 

ab 



We also invoke the following result Lemma 2, p354] (and see also p6|): 



{I (^sti) [{I ®q'''^) a] = KI, (101) 

for some constant K depending on the normalisation of a, where str denotes the 
supertrace. By expanding the LHS of ( p.01| ), we shall find an expression for 0~ 
sufficient to satisfy ( |100| ). In a sense, the result ( |101| ) amounts to satisfying the 



requirement that our invariant is an invariant of ambient isotopy. 



To summarise, we shall find: 



i^T = i-f (102) 
which, as = (13^) \ yields: 

so that we shall require the evaluation of the components 7r(g^^'''')^. 



To begin the expansion of the LHS of (101), we have: 



ijkl 



SO we must deduce the form of [(/ ® q '^^'') crYf,\- To this end, we have: 

\ im / \ jn 



i]mn 



i.e. in components: 

a ®?-'''')L'. = (103) 

Also: 

^ = E (^)rr (e'". ® e-,) , (104) 



mnkl 



TO 



Let A and B be general (ungraded) rank 4 tensors, viz A = J2ijmn ^In n i^^m ® e-'n) 
and B = Y.mnki (e™fc ® e",), we have: 

ij'fcZmn 

as, for eac/i m (i.e. without summation over m), we have e^mC^'k = Thus: 

Where Ais I ® q~^^p and 5 is a, we have the components: 

[(/ ® E^™(^""^)n(-)rz"- (106) 

mn 

If M is a rank 2 tensor (i.e. a matrix), viz: M = Mje*j, then we have: 

str(M) ^ E(-)'''^i' 
i 

hence if A is a rank 4 tensor, viz A = '^^jj^i A^f^\{e\ ® e^i): 
(/®str)(A) = E(-)^^''^l>V 

ijk 

i.e. componentwise: 

[(/®str)(A)]; = E(-)'^'^^U- (107) 



Applying the definition of the action of / ® str to the expression (^I (E) q '^^''^ a: 



b 

abm 



ab 



Now, by ( JlUlD , we have: {(/ ® str) [(/ ® q '^'^p) a] }^ = Kd^., for some constant K, 



that is we have: 

Y.(-f('i''"<(^)i: = ^^'- (108) 

ab 

Thus, we have K = 1. Choosing: 

iU-f'' = {-f\q-''^t (109) 
converts ( p.08| ) to (|100|) , and we have thus justified ( p.02| ). 



5.2.3 Evaluation of the Components (g '^^f) 



a 



Now we evaluate the components (g ^'*'')a- To this end, let H be the Cartan subal- 
gebra of gl{2\l), with dual the gl{2\l) root space H*. A basis for H* is: 

{£l=(l,0|0),£2 = (0,l|0),£3 = (0,0|l)}, 

on which we have an invariant bilinear form (-, ■) : H* x H* —>■ C defined by: 

{e,,e,) ^ (110) 

Now, gl{2\l) has the two simple positive roots ai = (1,-1 1 0) =61—62 (even) 
and ao = (0, 1 1 — 1) = 62 — 63 (odd), and another, nonsimple, odd positive root 
ai + ao = (1, I — 1) = 61 — £3. From these we have the half-sums of all even and 
all odd roots po = (f, — | I 0) and pi = (|, 1 1 — 1) respectively, hence their graded 
half-sum p = po — Pi is: 

p= (0,-l|l) = -£2 + ^3, 

from which we may proceed to determine the associated Cartan element hp. 



Firstly, as is a dual vector space to H*, we may implicitly define generators 
{-^%}?=i of H (i.e. gl{2\l) Cartan generators) by the action: 



E\{e,) ^ z, J = 1,2,3. (Ill) 

Next, we define hp ^ H to be such that: 

hp{ai) = {p,ai), i = 0,l. (112) 

Examination of the definition of p shows that for the two simple roots ao = 62 — £3 
and ai = £1 — £2, we have: 

{p,ai) = ^{ai,ai), « = 0, 1, 

thus ( |1 12|) may be more usefully written as: 

hp{ai) = ^{ai,ai), i = 0,l. (113) 

For gl{2\l) (as opposed to s/(2|l), for which it is unique), hp is in fact only determined 
up to addition of multiples of the first order Casimir C = E^i + + -£''^3, which 
satisfies C{ai) = for i = 0, 1. Thus, in setting hp = (3iE^i + /52-E'^2 + PsE^s, for 
complex coefficients Pi to be determined by ( |113| ); we may choose /?i = 0. Thus, 
using: 

(111 



we may write: 



hp{a^ = hp{e2) - hp{e?), =^2-^3 
hp{ai) = hp{ei) - hp{e2) = Pi - ^2- 



This information may be combined with data obtained by expansion of (TT^): 



hp{ao) = i(ao,ao) = 1(^2 - £3,^2 - £3) 1(1 - 1) = 
hp{ai) = i(ai,ai) = - £2,^1 - £2) ® 1(1 + 1) = 1 



to yield a linear system for the Pf. 



Thus, we have: P2 = = Pi — ^ = —1, and so, as /?i = 0: 



hp = —E^2 — E^2.- 

Note that this is expressly not + -£'^3 as might be naively assumed from 

inspection of p. 

The same hp holds for Uq[gl{2\l)], hence we have: 

H H 1 



and: 



7r 



l)(l|+g-^|2)(2| + |3)(3| + g-^|4)(4|]- 

|1)(1| + g2(a+l) |2)(2| + g2{«+l) |3)(3| + ^2{a+2) J^^^^jj 

|1)(1| + |2)(2| + |3)(3| + g2(a+i) |4^^4| ^ 



i.e., expressed clS cl matrix: 



^2a 



,2a 



2(a+l) 



2(q+1) 



p2g,2 



p2g,2 



(114) 



Note that this result disagrees slightly with that presented in the published work [jTSl , 
where the definitions of |2) and |3) are interchanged with respect to the above. This 
disagreement will not affect the evaluation of the invariant. 



Lastly, ( |114|) allows us to use ( |109|) to state 



p2q2 



and hence VL = (15 ) ^ is: 



p- 



-p- 



-p ^ 



p "^q 



6 Implementation in Mathematica 



6.1 Construction of the Braid Generator 



The program Uqgl21 .m (see § p.l|) consists of a suite of functions that evaluate the 
exphcit representation of ex, starting from the basis vectors {"^j) of V ® V. It also 
contains auxiliary functions required to implement noncommutativity and a few 
other functions to construct objects related to a. 

The code is self-explanatory. To use it, the program must first be imported into a 
Mathematica session via a command line such as Get ["Uqgl21 .m"] . After the 
functions are loaded, typing sigma will evaluate a, and typing ungradedsigma will 
evaluate a. Other functions behave similarly, and are named in an obvious fashion. 
Perusal of the code listing in § [B.1| will confirm what functions are available. Note 
that once a has been evaluated within a given session, Mathematica remembers 
its value, and when called upon again, it will not have to be recomputed. 

Of some interest is the code CheckRu that computes R{u). Whilst the output from 
that function is correct, considerable effort went into converting screen-based output 



to the human readable form of (|9l|) . 



6.2 Evaluation of the Invariant 



The program LinksGouldlnvariant .m (see § p.2| ) consists of a suite of functions that 
evaluate the Links-Gould invariant for various links, starting from a (computed by 
a function within Uqgl21.m) and a set of abstract tensors. It has been designed 
to compute the Links-Gould invariant and the (unnormalised) bracket polynomial 
(and thence the Jones polynomial). As it depends on Uqgl21.m, both that file and 
LinksGouldlnvariant .m must be imported. 

The command Linklnvariant [LinksGould, Link] returns the polynomial invari- 
ant corresponding to the link Link, which may be a pretzel Pretzel [pp, qq, rr] 
or any or all of the links in KnotList, which is: 

Unknot, Hopf Link, Trefoil, FigureEight, 
Cinquef oil, FiveTwo, WhiteheadLink, 
SixOne, SixTwo, SixThree, 
Septf oil, SevenTwo, 

EightSeventeen, NineFortyTwo, TenFortyEight, KT, KTI. 

The program also contains several auxiliary functions and a few other functions to 
evaluate the properties of the invariant. In particular, we have: 

• YangBaxterChecker [Invariant] asks if the a being used for the invariant 
satisfies the quantum Yang-Baxter equation. 

• LGCantDetectNoninvertibilityQ [Link] answers the question: "Is is that 
LG cannot detect the noninvertibility of Link (if it is noninvertible)?" by 
checking whether LG for Link displays the symmetry that denies it the ability 
to distinguish Link from its inverse. Link must be an element of KnotList, 
or be a Pretzel [pp, qq, rr] . 

• LGDetectsChiralityof LinkQ answers the question: "Does LG detect the chi- 
rality of Link (assuming it is chiral)?" by checking whether LG for Link 
is not palindromic. Again, Link must be an element of KnotList, or be a 
Pretzel [pp, qq, rr] . 

• Listof Pretzels [N] returns a list of the noninvertible pretzel knots {p,q,r) 
for all odd, distinct combinations of3^p<g<r^A^ (odd). 

• LGCantDetectNoninvertibilityof PretzelsQ [N] returns a list of answers to 
LGCantDetectNoninvertibilityQ [Pretzel [pp,qq, rr] ] for all the pretzels 
Pretzel [pp,qq,rr] in Listof Pretzels [N] . 

• LGDetectsChiralityof PretelsQ [N] returns a list of answers to the ques- 
tion LGDetectsChiralityof LinkQ [Pretzel [pp,qq,rr] ] , for all the pretzels 
Pretzel [pp,qq,rr] in Listof Pretzels [N] . 

• JonesPolynomial [Link] evaluates the Jones polynomial for the link in ques- 
tion. 

TO 



7 Experimental Results and Analysis 



7.1 Analysis of the Behaviour of the Invariant 

Before proceeding to evaluate the invariant for our test examples, we digress to 
consider its properties so as to be able to make predictions about the evaluations. We 
are interested in whether LG distinguishes the chirality of links, the noninvertibility 
of (single-component) knots and mutations of knots. 



7.1.1 Chirality: The Effect of Reflection on LGk 

Say that we have a link K with associated abstract tensor and invariant LGk- If 
we replace K with its reflection K* , every positive (respectively negative) crossing in 
K will have been replaced by the equivalent negative (respectively positive) crossing 
in K* , thus, T^* will be the same as except for the interchange a ^ a~^. The 
caps and cups will remain unchanged, as will the directions of the arrows 
(recall that chirality is independent of orientation). 



From the uniqueness of the universal R matrix for any quantum (super) algebra [33 
the following relation holds (for appropriate normalisation): 

R-\q)=Riq-'), 



hence, as a = PR, we have: 



a-\q) = Pa{q-')R 

(In the case of the Links-Gould invariant, R of course has the extra variable a, but 
this doesn't affect the above.) Thus, up to a basis transformation, a and are 
interchangeable by the change of variable q ^ q"^ . It then follows that the invariant 
for K* is obtainable from that of K by the same change of variable, thus we have 
proven: 



Proposition 1 

If K is amphichiral then the polynomial LGxiq) is palindromic, i.e. invariant 
under the mapping q ^ q~^ (and hence p p~^). 



From the polynomial for a link we may immediately write down the polynomial 
for K*: 

LGK4q,p) = LGK{q-\p-') (115) 
As we have: 

K = K* LGK{q,p)=LGK*{q,p), 
then we have, conversely, that: 

LGK{q,p)^LGK'-{q,p) ^ K ^ K\ (116) 

i.e. if the polynomials corresponding to K and K* are distinct, then K must be 
chiral. Using the identity (|115|) , the test of (|116|) becomes: 

LGK{q,p)^LGK{q-\p-') ^ K ^ K\ 
hence we have demonstrated: 

Proposition 2 

If LGxiq-iP) is not palindromic, then LG detects the chirality of K. 



en 



7.1.2 The Effect of Inversion on LGk 



Replacing K witli amounts to reversing every arrow in K. T^-i will thus be the 
same as Tk except for the following changes, best expressed in terms of the auxiliary 
tensors. For the crossings, where X is either a or cr"^, interchange X with X^ and 
Xi with Xr'i and for the caps and cups, interchange only the signs, i.e. interchange 
fi± with and 0± with O^. 

This has the effect that Tk is replaced by T^-i, the dual tensor acting on the 
dual space Recalling that the tensors representing (1,1) tangles act as scalar 



multiples of the identity on V , the dual tensor has exactly the same form, thus: 



Proposition 3 

A knot invariant derived from an irreducible representation of a quantum 
(super) algebra is unable to detect inversion. 

Now observe that the representation of [/q[5f/(2|l)] acting on the dual module V* 
is given by the replacement a ^ —{a + 1) (with an appropriate redefinition of the 
Cartan elements). Thus for a given (1,1) tangle with invariant LGk{(1iP)i the 
invariant LGj^-i of its inverse is obtained as: 

LGK~iiq,p) = LGKiq,q'^p''^)- 

However, in view of Proposition such an invariant is unable to detect inversion, 
hence we must have LGk-^ = LGk, from which it follows that: 



Proposition 4 

LGxiqiP) enjoys the symmetry property: 

LGK{q,p) = LGk {q,q-^p-') • (117) 



7.1.3 The Effect of Mutation on LGk 



The question of whether our invariant is able to distinguish mutants is answerable 
in the negative. Theorem 5 of states that (for quantum algebras) if the modules 
occurring in the decomposition oi V each have unit multiplicity, as indeed (|64D 



shows in our case, then the invariant is unable to detect mutations. Whilst this was 



proved in |^ for the case of quantum algebras, the extension to the case of quantum 
superalgebras is quite straightforward. We make a proposition out of it: 



Proposition 5 

The Links-Gould invariant is unable to distinguish between mutants. 



Q1 



7.2 Invariant Evaluations for the Ordinary Links 



Evaluations of LGxiq, P) collection of ordinary links are presented in Table 0. 

We have, as intended, that LGqi = 1- Observe that LGq.^ are palindromic 

as 4i and 63 are amphichiral, whilst LGk is not palindromic for any other links as 
the rest are chiral. Observe also that the invariant is actually a polynomial in (and 
hence p^) rather than just q (and p), thus we might write LGk{Q = (f^P — P^) 
rather than LGk^I^p)- 



K 


LGK{,q,p) 


2? 


—1 + — + p^g^ 


3i 


1 + p~^ — + — p"^g^ — p^g^ — p^q^ + p'^g'' 


4i 


7 + (p-4g-2 + /^2) _ 3 (p-2 + _ 3 (p-2^-2 ^ p2g2) ^ 2 (g-2 + g2) 


5i 


1 _l_ p-8 _ p-6 _^ _ p-2 _^ 2g2 _ ^-6^2 _^ 2p~^g2 - 2p^^g2 - p^q^ + 2g^ 
_p-2qi — 2p^g^ + p^g^ — p^g^ + 2p^g^ — p^g^ — p^g^ + p^g*^ 


52 


3 + 3p-^ - 5p2 + 10g2 + p-4g2 - 6p-2g2 _ 5p2^2 ^ 4^4 _ ^-2^4 _ gp2^4 

A A (-1 A f\ 

+3p g — p g + p g 


5? 


-10 + p-6g-2 - 3p-^ - 3p-^g-2 + 4p-2g-2 + 9^-2 _ 2g-2 

-8g2 + 2p^2^2 ^ gp2^2 ^ 4^2 ^ 2p2g4 _ 3p4^2 _ 3p4^4 ^ ^6^4 


61 


1 7 _L 4„— 4 _|_ Q„— 4„— 2 0^— 2„— 4 1 r\^—2„—2 7r,^2 i o„— 4 i 1 fi^— 2 

i ( + p g + op g — op g — iup g — / p + zq -\- rug 
+4g2 - 3p^g-2 - lOp^ - 7p^q^ + p"^ + 3p'^g2 


62 


11 - 3p-^ + 9p-^ - 12p-2 - 4p2 + 2g-2 + p-8g-2 - Sp'^g-^ + 4p-^g-2 
-4p-2g-2 + 14g2 + 2p-^g2 - 8p-2g2 - 12p2g2 + 4p^g2 + 2g^ - Sp^g^ 
+9p^g^ - 3p6g^ + 2p^g6 - Sp'^g^ + p^g^ 


63 


25 + {p-^q-^ + p^q^) + 10 (g-^ + g2) + H (^-4^-2 ^ ^4^2) ^ 4 (^-4 ^ ^4^ 

+4 (p-4g-4 + p4g4) - 2 {p-^q^ + p2g-2) _ 16 (p-2 + p^) 

-16 (p-2g-2 + p2g2) _ 2 (p-2g-4 + ^2^4) _ 3 (^-6^-2 ^ ^6^2) 

-3 (p-6g~^ + p6g^) 


7i 


1 + p"^^ - p"^° + p^^ - p~^ + p~^ - p~^ + 2g2 - p^°g2 + 2p-'^g2 - 2p-''g2 

+2p-4g2 - 2p-2g2 _ p2^2 ^ 2g4 _ p-6g4 ^ 2p-4g4 - 2p-2g4 _ 2p2g4 

+p^g^ + 2g6 - p-2g6 _ 2p2g6 + 2p^g6 - p^g^ - p^g^ + 2p^g8 - 2p6g* 

_j_p8g8 pGglO _j_ 2,p^q^^ plOglO pl0gl2 _^ ^12^12 


72 


5 + 5p~^ - 9p-2 + 20g2 + 3p-^g2 - 14p-2g2 _ 9p2^2 ^ ^4^4 ^ ^-4^4 
—Qp^'^q'^ — Mp^g"^ + 5p^g^ + 4g^ — p"^g^ — 6p^g^ + 3p'^g^ — p^g^ + p'^g® 



Table 2: LGK{q,p), evaluated for various oriented links K. Data for 817 are presented 
in §|773| and for 942 and IO48 in WA- 



Evaluations of the Jones polynomials Vxit) for our example links are presented in 
Table 1^. These polynomials agree with published data Their correctness 

confirms that the evaluations of the Links-Gould invariant are correct. 



K 




1 


1 5 


3i 


t + t^-t^ 


4i 




5i 




52 


t-t^ + 2t^ -t* + - 


1 


-f-l + - 2U + tt - 2tt + ti 


6i 


t-* - t-'^ + _ 2t-i + 2 - t + t2 


62 


t-^ -l + 2t- 2t^ + - + 


63 


- + t3) + 2 + _ 2 + t) + 3 


7i 




72 


t - + 2^3 - + 2t^ - + t7 _ ^8 




017 


[t + t ) — 6 [t + t j + 5 (r +t) — b[t + rj + 7 


942 




IO48 


- (t-^ + t^) + 2 + t^) - 4 + t3) + 6 (r2 + t2) 

-7(ri + t) + 9 


KT, KTI 


t-6 _ 2t-5 + 2t'4 - 2t-^ + + 2t - 2t2 + 2t^ - 



Table 3: Vxit), evaluated for various oriented links K. 



7.3 The Noninvertibility of 817 is not Detected 



Recall that 817 is the smallest noninvertible knot. We find: 

LG8,,{q,p) = 

139 + {p-'^q-^ + p^^q^) - 4 {p-'^q-' + p'^q') - 4 {p~''q~' + p''q') 
+7 {p~^q~^ + p^q^) + 18 [p'^q-^ + p^q'^) + 7 {p~^q~^ + p^q^) 
-7 [p-^q-^ + /g^) - 36 [p-^'^ + /g^) - 36 {p'^q''^ + /g^) 
-7 (p-' + p') + 3 + p^g6) + 40 (p'^g"^ + p^g^) 

+82 {p'^q-^ + p'g2) + 40 [p-^ + p') + 3 (p-^g' + p^g"') 
-22 (p'^g"^ + p^g^) - 102 (p"^g"^ + p^g^) - 102 [p'^ + p^) 
-22 (p-2g2 + p2g-2) + 4 (g-^ + g^) + 68 (g"' + g^') . 

We observe the invariance of Proposition ^, viz LG^^^i^q.p) = LGs^^iq, q^^p^^), 
illustrating that LG doesn't detect the noninvertibility of 817. Also note that as 817 
is amphichiral, the polynomial is palindromic, as predicted by Proposition |1]. 

7.4 The ChiraUty of 942 and IO48 is Detected 
We have: 

LGg,,{q,p) = 

3 + p-^q-^ - 2p-\-^ - 2p-^q-^ + p'^g-^ + 3p-^g-^ + p~^g~^ + p"^ 

-2-4 -2-2 o -2 o -2 2 I p 2 I r) 4 2 -2 2 o 2 2 

— p g — p g — 3p — 3p q + bq + 2q — p q — p — 6p q 
~3p'^q'^ + p^g-2 + 3p^ + p^g^ + p^g^ - 2p'' - 2p^q'^ + p^q^ 
LGio^,{q,p) = 

165 + 5p"^ - 25p~^ + 68p"^ - 129p"^ - 132p2 + 67p^ - 22p^ + 4p^ 

+p-16g-8 _ 3p-14g-8 ^ 4p-12g-8 _ 4p-10^-8 ^ 4^-8^-8 _ 3^-6^-8 

-3p-^^g"^ + 12p-^\-^ - 21p~i°g~^ + 24p-®g-'' - 22p-6g-^ + 13p"^g-^ 
-3p"2g"^ + 16g"^ + 5p"^^g"^ - 23p~^°g"^ + 50p"*g~^ - 69p"^g"^ 
+67p"^g~^ - 43p~2g"^ - 3p2g-^ + 94g~2 _ 6p~^^q-^ + 29p-*g-2 
-72p-^g-2 + 119p-^g"2 - 132p~2g"2 _ 43^2^-2 + 13^4^-2 _ 2/^-2 

+88g2 - 2p"^g2 + 12p-^g2 - 39p"^g^ - 129p2g^ + 119p^g2 - 69p^g2 
+24pSg2 - 4p^°g2 + 12g^ - 2p~^q'^ - 39p2g^ + 68p^g^ - 72p6g^ + 50p®g^ 
-21p^°g^ + 4p^2g^ - 2p2g^ + 12p^g^ - 25p^g^ + 29p^g^ - 23p^°g^ 
+12p^^q^ - 3p^^g^ - 2p^g^ + 5p*g^ - 6p^°g^ + Sp^^g* - 3p^^g* +p^^g^ 



As neither of these polynomials are palindromic, LG distinguishes the chirality of 
both these knots. 

Q A 



7.5 The Noninvertible Pretzels are not Distinguished 



Experiments show that the Links-Gould invariant for this class of noninvertible 
knots always displays the symmetry of ( |117| ), for all p,q,r ^ 67. This amounts to 
5456 knots0 (sheer bloody-mindedness!), the smallest being the (3,5,7) pretzel, a 
knot of 3 + 5 + 7 = 15 crossings, and the largest being the (63, 65, 67) pretzel, a knot 
of 63 + 65 + 67 = 195 crossings. This illustrates the assertion of Proposition ^ that 
our invariant cannot detect the noninvertibility of any knot. The smallest of these 
polynomials is: 



57 + 57p'^ - 113p-2 + 308g^ + SAq^p'^ - 279q'^p-^ - llSp'^q^ + 468g^ 
+83g^p-^ - 329q^p-^ - 279p^q^ + hlp^q^ + 464g^ + 60g^p-^ - 279q^p-^ 
-32Vg6 + 8ApY + 338g^ + 30q^p-^ - 172q^p-^ - 279p^q^ + 83p^g^ 
+ 174g^° + 9g^V"^ - 71g^V"^ - I72p^q^^ + 60/g^° + 56^^^ + q^'^p-^ 
-Wp-' - 71pV' + 30p'q'' + 8q'' - q'^p^' - IQp^^ + 9pW' 
-p'q'^ + pY''. 



We also discover that all these 5456 pretzels are chiral as none are palindromic, 
lending weight to Trotter's assertion that all such pretzels are chiral. 



7.6 The Kinoshita— Terasaka Mutants are not Distinguished 

The polynomials for both mutants are: 

LGKT{q,p) = 

-23 - p-^q~^ - p-^q-^ + 2p-^q-^ + p-^q-^ - p-^ + p-^q-^ + Gp'^q'^ 

-?,p~^q'^ - 9pS~^ + 2p-^ + 3p-^g2 - 7p~^q-^ - 7p~^q^^ + I8p'^q~^ 

+9p-^ - llp-'^q^ - 2p~^q^ + 2q-^ + Uq-^ - Sq'^ + 6q^ + lOq^ - 7p^q-^ 
-7p^q-^ + 18p2 + 9/g2 _ 11^2^4 _ 2p2^6 ^ ^4^-4 ^ Q^i^~2 _ 3^4 _ ^^^^2 

+2/g^ + 3p^g^ - p'^q''^ - / + 2p'^q^ + p^g^ - p^g^ 



hence LG does not distinguish between these mutants. As predicted by the theorem 
of (i.e. Proposition ^), the tensors KTA and KTA' are in fact identical, which 



explains why the pair of mutants yield the same invariant. 



^ More generally, if we sample for all p,q,r ^ N, for some odd TV ^ 7, we will have tested 
{N — 1)(-/V — 3)(iV — 5)/48 pretzels, the largest of which is a knot of 3N — 6 crossings. 



8 Conclusions and Extensions 



8.1 Conclusions 



The Links-Gould invariant is at least as powerful as the HOMFLY and Kauffman 
(two- variable) invariants in the following sense. Neither of those invariants can detect 
the chirality of 942, although they can do so for all preceding knots. Beyond 942, 
there are other knots for which they fail to detect chirality. IO48 is the next example 
for which the HOMFLY invariant fails to detect chirality (although the Kauffman 
does detect it). However, although our invariant does detect the chirality of 942 
and IO48, there is as yet no evidence to tell us whether it detects the chirality of all 
preceding knots, that is, for all we know, it may fail on 941 (which is chiral, according 
to p385], p2| , p265]). This information is summarised in Table ^. 



Knot 


HOMFLY 


Kauffman 


Links-Gould 


3i to 72 


Yes 


Yes 


Yes 


73 to 941 


Yes 


Yes 


? 


942 


No 


No 


Yes 


943 to IO47 


Yes 


Yes 


? 


1048 


No 


Yes 


Yes 


IO49 and above 


sometimes 


sometimes 


? 



Table 4: Chirality detection properties of two-variable link invariants. 

We have as yet found no chiral knot for which LG fails to detect its chirality, so 
whether our invariant is a complete invariant for chirality is an open question. 



Recall that the method used to evaluate our invariant for a particular link involves 
the deduction of an appropriate abstract tensor that is suitable for use with efficient 
shortcuts using auxiliary tensors, and this requires the labour-intensive drawing of 
a suitable presentation (i.e. finding an appropriate quasi- Morse function). It is thus 
infeasible to use this method to systematically evaluate the invariant for all links up 
to any particular point. This situation is in marked contrast to that for the HOMFLY 
and Kauffman polynomials, for which extensive tables have been (automatically) 
computed. The question of whether the Links-Gould invariant is superior to the 
other two- variable invariants is therefore currently unanswerable by experiment. To 
this end, it would be useful to conduct research to develop a method to automate the 
deduction of the abstract tensor directly from a braid presentation. (The abstract 
tensor would then be applicable to the evaluation of any other invariants based on 



state models.) Work on this program is in preparation ||TT 



As explained in Proposition ^ no invariant of this typeQ can be expected to detect 
inversion. Our early hopes that our invariant might do so were shattered by exper- 
imental results; indeed this led us to the literature on invariants constructed from 
quantum algebras, which led us to propose that the same considerations naturally 
carry over to quantum superalgebras. 

Also, as explained in Proposition |], no state space invariant of this type^ can be 
expected to distinguish mutants; again negative experimental results caused us to 
make the considerations leading to the theoretical proof of this statement. 



8.2 Extension — A ?7q[5'/(3|l)] Two- Variable Link Invariant? 



Having investigated the Links-Gould invariant based on the f/q[5fZ(2|l)] family of 
representations labelled (0,0|a), it is immediately attractive to ask if other two- 
variable invariants may be constructed from the families (Om | a) of Ug[gl{m\l)] 
representations, for general m ^ 2. The immediate answer is that yes, of course 
they exist. Their actual evaluation would require a lot more computation as the R 
matrices would be larger, but the process would in principle be the same. In practice, 
explicit construction of the R matrices is hampered by the algebraic complexities 
involved in finding an explicit orthogonal basis for the tensor product module. 



8.2.1 An Orthonormal Basis for V 



We illustrate the situation with the example of the ^^[(^/(Sll)] family of represen- 
tations (0,0,0 I a). This family is 8 dimensional (cf. the family for [/g[5f/(2|l)] is 4 
dimensional), which means that the R matrix has 8^ = 4096 rather than 4^ = 256 
(albeit mostly zero) entries, so the computational effort in evaluating abstract ten- 
sors would be orders of magnitude larger. 



A basis for the decomposition of is provided in |T^, but this is not orthogonal. 

The decomposition involves two modules of dimension 8, for which a fully orthogonal 
basis is given, and two modules of dimension 24, of which a basis for only one (called 
V2) is provided, and that basis is not orthogonal. Using the Gram-Schmidt process 
to (manually!) orthogonalise that basis for V2 yields reasonably tractable expressions 
for all but the last vector (originally called |^'24))- Attempting to normalise this basis 
would further complicate the expressions involved, and it appears that the resulting 
R matrix would contain expressions similar to Y (see (PP|)) in our R matrix, although 
somewhat more complicated. The complexity of Y doesn't propagate through to 
evaluations of LG, and we anticipate that evaluations of the Uq[gl{3\l) invariant 
would behave similarly. 



^ That is, an invariant based on a state space model where the crossing tensor a is based on 
the R matrix associated with a quantum (super)algebra representation. 

^ That is, an invariant based on a state space model where the crossing tensor a is based on the 
R matrix associated with a quantum (super)algebra representation in which the decomposition of 
the tensor product contains only unit multiplicities. 



QO 



8.2.2 Sour Grapes 



Perhaps there is nothing lost by the lack of success in constructing a state model 
for an invariant based on [/g[^/(3|l)] family of representations (0,0,0 | a). Such an 
invariant would have the same inability to detect inversion and mutation as does 
LG, and there is no particular reason to expect that it would always detect chirality. 



8.3 Automatic Construction of Uqlgl^mln)] Representations 

The process used to build the explicit Uq[gl{2\l)] representation was presented in 
rigorous detail so that it might be analysed in order to consider the automatic con- 
struction of Uq[gl{m\n)] (generally non-parametric) representations. (A description 
of Uq[gl{m\n)] is provided in §0.) 

A general element of f/g[5'Z(m|n)] is a sum of (complex) multiples of products of 
Uqlgllmln)] generators. As Uqlgllmln)] is a (graded) vector space, we may meaning- 
fully call these general elements "vectors" . 

So, given a vector involving weighted strings of Uq[gl{m\n)] generators, where each 
string may contain nonsimple generators, how do we convert this to some kind of 
normal ordering? Well, firstly we must (recursively) replace every nonsimple gen- 
erator in the vector with a combination of simple generators using ( |119D , until our 
vector has only simple generators. Secondly, we must (recursively) examine every 
product of two generators "Ai?" in each summand in the vector. If the pair of gener- 



ators are disordered (i.e. A > B),we use the relations of § |A.3| to replace the product 
with another, ordered string. 

This process provides us with all the necessary machinery for the manipulation of 
vectors encountered when constructing concrete representations of Uq[gl{m\n)]. 



on 



A The Quantum Superalgebra Uq[gl{'m\n)] 



The origin of the quantum superalgebra Uq[gl{m\n)] (we intend 171,11 E Z+) as a 
quantum deformation of gl{m\n) is described in [^,0,^, 
by Chari and Pressley p], see §6.5]. Here, we set the context for § |4.2| (with obvious 
extensions of the definitions of the graded commutator and q bracket). The mate- 
rial is taken from |5^, ppl237-1238], except that we have substituted the notation 
^(-j'^i^a modified the definition of the coproduct and antipode. 



50,51], and in the book 



A.l Z2 Grading on Uq[gl{'m\'n)] Indices 

Firstly, we define a Z2 grading [■] : {1, . . . ,m + n} ^ Z2 on the m + n gl {m\n) 
indices: 

|. 1 _ J if a = 1, . . . , m (even) 
[1 if a = m + 1, . . . ,171 + n (odd). 

Throughout, we shall use dummy indices a,b = 1, . . . ,in + n where meaningful. 



A. 2 Uq[gl{m\n)] Generators 
A set of generators for Uq[gl{m\n)] is: 

's^a^ a = 1, . . . ,m + n (m + n Cartan) 

E°'iy, a < b = 1, . . . ,m + n ((m + n)(m + n — l)/2 raising) 

E^a, a < b = 1, . . . ,m + n ((m + n)(m + n — l)/2 lowering). 

We write the Cartan generators as g'-"-'' they might also be written as g^""; 
the notation including the parity factor ensures consistency with the literature and 
what follows below. More significantly, we are not writing the Cartan generators as 
E^'a (as per gl{m\n)). 

On all the (m + n)'^ generators we define a natural Z2 grading in terms of the grading 
on the indices: 

[E\] = [a] + [b] (mod 2), 

where the definition of the exponential yields [q^""-] = [E°'a\ = 0. The product of 
homogeneous x,y E Ug[gl{m\n)] has degree: 

[xy] = [x] + [y] (mod 2) . 



A. 2.1 Uq[gl{m\n)] Simple Generators 

The Uglgllmln)] simple generators are a subset of 3{m + n) — 2 elements of the above 
generators: 

'-B^a^ a = 1, . . . ,m + n (m + n Cartan) 

E"'a+i, a,a + 1 = 1, . . . ,m + n (m + n — 1 raising) 

E°'~^^a, a,a + 1 = 1, . . . ,m + n (m + n — 1 lowering). 

where by the expression E°-a+i, a,a + 1 = 1, . . . , m + n, we of course intend to 
exclude the nonmeaningful E'^i and £''"^"m+„,+i. Note that there are only two odd 
simple generators: E"^rn+i (raising) and E"^~^^m (lowering). The fact that there are 
m + n — 1 simple raising generators indicates that Uq[gl{m\n)] has rank m + n — 1. 



A. 2. 2 Uq[gl{m\n)] Nonsimple Generators 

In the 5f/(m|n) case, the remaining nonsimple (non Cartan) generators satisfy the 
same commutation relations as do the simple generators. Writing e°'b as the (m + n 



,2 



gl(m\n) generators, to parallel the E^-i, of Uq[gl{m\n)], we have: 

[e\ e'^,] = - (-)(H+M)(W+H)5a^e^,, (118) 

and this is true for all generators e"?,, e'^d, not just the simple ones. To be sure, Zhang 
pl238] obtains these generators by taking the limit as g ^ 1 of Uq[gl{m\n)] to 
yield f/[(7/(m|n)], with simple generators: 

e'^a+i = limq^i E°'a+i, a = 1, . . . ,m + n — 1 (m + n — 1 raising) 
e"a_i = limq^i E"'a-i, a = 2, . . . ,m + n (m + n — 1 lowering) 

e°'a — liHig_>i [E"'a]q, tt = 1, . . . ,m + u {m + n Cartan). 

(Of course, we have lim^^i g^^"" = 1.) Having defined these, the gl{m\n) commu- 
tation relations ( [118| ) serve to define the nonsimple generators. 

The situation is different in Uq[gl{m\n)]. The remaining (nonsimple) generators do 
not satisfy the same commutation relations, in fact they are recursively defined in 
terms of the simple generators ||5^, pl971, (3)] and [^, pl238, (2)]: 



(a) E^b = E''cE%- q-^-'^^'^E%E''c (raising) 

(b) E\ ^ EKE^-q^'^^'^E^E^ (lowering) 



(119) 



where a < c < 6. To evaluate a graded commutator involving a nonsimple generator, 
that generator must first be expanded using (|119|) , and then the graded commutator 
expanded by linearity. 



no 



A. 3 Uq[gl{'m\n)] Relations 

For the simple Uglgllmln)] generators, we have the following Uq[gl{m\n)] relations: 

• The Cartan generators all commute: 

I T7a I Tpb I 776 I jpa jpa rpa 

n "-a * = a " n n~ " = 1 

• The Cartan generators commute with the simple raising and lowering genera- 
tors in the following manner: 

^(_)H^.^^,^^^^_(_)W^.„ = (120) 
From ( |12CI| ), we have the following useful interchange: 

Replacing q with in ( |121| ) yields the equivalent: 

q ( ) ^ b±i = q ' ^>'±^)E\±iq ^ > ^ ^ 

• The squares of the odd simple generators are zero: 

It may be shown that this implies that the squares of the nonsimple odd 
generators are also zero. 

• The non-Cartan generators satisfy the following commutation relations (this 
is the really interesting part!): 



+1 



We also have, for |a — 6| > 1, the commutations: 

-C/ a+l-C/ b+1 — ^ h+li^ a+l ana £/ 6 — -C/ q. 



no 



There are two additional relations, known as the Serre relations. Their inclu- 
sion ensures that the algebra is reduced enough to be simple. For a ^ m: 

{E'^ a+lf E''^^ a±l+l — (Q' + q'^) E"- a+lE""^^ a±l+lE°-' a+1 

+E°-^^ a±l+l{E'' a+l) = 

+E-^'+\^^{E-+\f = 0. 
To complement these, for a limited set of odd generators, we have: 

[E"^~^m+2, E"^rn+l] = [E"^^"^ m-l, E"^'^^ rn] =0, (122) 

where, as all the elements are odd, the graded commutator devolves to a simple 
anticommutator, viz (|122|) may be written as the exchanges: 



rpm—l Z7"TS TTim rpm—1 

^ m+1 — m+l-C/ m+2i 

7-im+2 7-im+l 7-im+l rp'm+2 



m—1- 



Note the presence of the nonsimple generators E"^ ^m+2 and E'^'^'^m-i here. 
The index range ensures that these relations are irrelevant to f/q[(7/(2|l)], which 



is why we didn't encounter them in §4.2.3 



It must be emphasised that expressions involving nonsimple generators may in gen- 
eral only be manipulated by prior expansion using (|119|) . 



A. 4 Uq[gl{'m\n)] as a Hopf Superalgebra 

Ug [gl{m\n)] may be regarded as a Hopf superalgebra when equipped with the fol- 
lowing coproduct A, counit e and antipode S structures [|^, pl238]. It is in fact 
quasitriangular (i.e. it possesses an R matrix). 



A. 4.1 Coproduct A 

We define a coproduct (a.k.a. comultiplication) structure, which is a Z2 graded 
algebra homomorphism A : Ug [gl{m\n)] Uq [gl{m\ny\ ® Uq [gl{m\ny\ by: 

(c) A(g±(-)""^"'') = ® g±(-)'"'^'"^ 

and extended to an algebra homomorphism on Uq [5f/(m|n)]. A being graded means 
that it preserves grading, viz for homogeneous x & Uq [gl{m\n)], we have [A(x)] = [x]. 
A being a homomorphism means that A(l) = 1 1, and for x,y & Uq [gl{m\n)]: 

A{xy) = A{x)A{y). 



As well as this standard coproduct, there exists another possible coproduct structure: 
A, defined by A = T ■ A, where the twist map T, an operator on the tensor product 
Uq[gl{m\n)] Uq[gl{m\n)], is defined for homogeneous x,y E Uq[gl{m\n)] by: 



A. 4.2 Counit e 

We define a counit e : Uq [gl{m\n)] C, also a Z2 graded algebra homomorphism: 
and extended to an algebra homomorphism on all of Uq [gl{m\n)]. 



A. 4. 3 Antipode S 



We define an antipode S : Uq [gl{m\n)] — > Uq[gl{m\n)], which is a Z2 graded algebra 
an^i homomorphism (actually an anti automorphism): 

and extended as an algebra homomorphism on all of Uq [gl{m\n)]. Because S is an 
anii automorphism, i.e. for homogeneous x,y E Uq[gl{m\n)]: 

S{xy) = {-p^'^S{y)S{x), 
we may deduce that 5* (1) = 1. 



B Mathematica Code 



B.l Uqgl21.m 



(* 



Uqgl21.m 



Mathematica code to build projectors and the braid generator. 

We are interested in the 4 dimensional representation LA = (0,0 1 a) 

of the Quantum Superalgebra U_q[gl(2|l)] . 

David De Wit 



(* Basis vectors of V \otimes V *) 

ket[l, 1] = Kron[ket[l], ket[l]]; 

ket[l, 2] = 

(q-a + q-(-a))-(-l/2) * 

(q-(a/2) * Kron[ket[l], ket[2]] + q-(-a/2) * Kron[ket[2], ket[l]]); 

ket[l, 3] = 

(q-a + q-(-a))-(-l/2) * 

(q-(a/2) * Kron[ket[l], ket[3]] + q-(-a/2) * Kron[ket[3], ket[l]]); 

ket[l, 4] = 

(q-a + q-(-a))-(-l/2) * 
qbracketCq, 2 a + l]-(-l/2) * 
( 

Qbracket[q, a + l]-(l/2) * 

(q-a * Kron[ket[l], ket[4]] + q-(-a) * Kron[ket[4], ket[l]]) 
Qbracket[q, a] -(1/2) * 

(q-(l/2) * Kron[ket[3], ket[2]] - q-(-l/2) * Kron[ket[2], ket[3]]) 

); 



20 
20 
08 



May 1996 
March 1998 
October 1998 



3 June 1996 

27 March 1998 
20 November 1998 



*) 



ket[3, 1] = 

(q-(a+l) + q-(-a-l))-(-l/2) * 
QbracketCq, 2 a + l]'(-l/2) * 
( 

Qbracket[q, a] "(1/2) * 

(q'(a+l) * Kron[ket[4], ket[l]] + q'(-a-l) * Kron[ket[l], ket[4]]) 
+ 

Qbracket[q, a + l]"(l/2) * 

(q-(l/2) * Kron[ket[3] , ket[2]] - q-(-l/2) * Kron[ket[2], ket [3] ] ) 

); 

ket [3, 2] = 

(q-(a+l) + q-(-a-l))-(-l/2) * 
( 

q"((a+l)/2) * Kron[ket[4], ket [2] ] 
+ 

q-(-(a+l)/2) * Kron [ket [2] , ket [4] ] 

); 

ket [3, 3] = 

(q-(a+l) + q-(-a-l))-(-l/2) * 
( 

q-((a+l)/2) * Kron[ket[4], ket[3]] 
+ 

q"(-(a+l)/2) * Kron [ket [3], ket [4]] 

); 

ket [3, 4] = Kron [ket [4], ket [4]]; 

Grading [ket [k_] ] := (Grading [ket [k] ] = If [MemberQ [{1 ,4} ,k] ,0, 1] ) ; 

Grading [bra [k_] ] := (Grading [bra [k]] = Grading [ket [k] ]) ; 

DimTPSubModule[l] := 4; 

DimTPSubModule [3] := 4; 

bra[k_, j_] := 
(bra[k, j] = 
Expand [ 

ket[k, j] //. 

Kron [ket [h_] , ket [i_]] :> 

(-1)" (Grading [ket [h]]*Grading [ket [i]]) * Kron [bra [h], bra[i]] 

] 

) /; 

(MemberQ [{1, 3}, k] && (j <= DimTPSubModule [k] )) ; 
(* *) 



no 



(* Utility functions *) 

KroneckerDelta[i_, j_] := If[i == j, 1, 0]; 

OperatorQ [x_NonCoinmutativeMultiply] := Map [OperatorQ , Apply [Or, x]]; 
OperatorQ [x_Times] := Map [OperatorQ , Apply [Or, x]]; 
OperatorQ [x_Plus] := Map [OperatorQ, Apply [Or, x]]; 

OperatorQ [ket [_] ] := True; OperatorQ [bra [_] ] := True; 

OperatorQ [Kron [_ , _] ] := True; 
OperatorQ [_] := False; 

SetAttributes [NonCommutativeMultiply, {Listable, Flat, Oneldentity}] ; 

Unprotect [NonCommutativeMultiply] ; 

NonCommutativeMultiply [Kron [A_, B_] , Kron[C_, D_]] := 
(-1)" (Grading [B] * Grading[C]) * Kron[A ** C, B ** D] ; 

(X_) ** ((alpha_) (Y_)) := alpha X ** Y /; NumberQ [alpha] ; 

(X_) ** (alpha_) := alpha X /; NumberQ [alpha] ; 

** X_ := /; OperatorQ [X] ; 

((a_)*(b_)) ** (c_) := a*b ** c /; ! OperatorQ [a] ; 

(c_) ** ((a_)*(b_)) := a*c ** b /; ! OperatorQ [a] ; 

((a_) + (b_)) ** (c_) := a ** c + b ** c; 

(a_) ** ((b_) + (c_)) := a ** b + a ** c; 

Protect [NonCommutativeMultiply] ; 

(* 

(* Build projectors: *) 
M = 4; 

Iden := Outer [Times, IdentityMatrix [M] , IdentityMatrix [M] ] ; 

Projector [K_] := 
(Projector [K] = 
Table [ 

Simplify [ 

Coefficient [ 

Sum[ket[K, j] ** bra[K, j] , {j ,DimTPSubModule [K] }] , 
Kron[ket[i] ** bra[j] , ket [k] ** bra[l]] 

] 

], 

{i,M},{j,M},{k,M>,{l,M} 
]) /; MemberQ[{l,3},K] 

Projector[2] := Simplify [Iden - (Projector [1] + Projector [3] )] ; 

(* 



nn 



(* Build the braid generator sigma and associated objects *) 



sigma := sigma = 
Module [ 
{ 

temp = 

q~(-2a) Projector[l] - Projector[2] + q"(2a+2) Projector [3] , 
ruleqX = {-1 - q~2 + q~(-2a) + q~(2+2a) :> X}, 
ruleqatop = 

{ 

q-(-a) :> p-(-l), q-(-2a) :>p-(-2), 

q'(l+a) :> q p, q'(2+2a) :> q~2 p~2 

} 

}, 

temp = temp //. Qbracket[q_, x_] :> (q~x - q~(-x))/(q - q~(-l)); 
temp = Simplify [Expand [Simplify [temp] ] //. ruleqX] ; 
temp = Simplify [PowerExpand [temp] ] ; 

temp [[4, 3, 1,2]] = -q Sqrt [Expand [(temp [[4, 3,1,2] ]/q)~2] //. ruleqX]; 
temp [[3, 4, 2,1]] = q Sqrt [Expand [(temp [ [3,4,2, 1] ] /q) ~2] //. ruleqX]; 
temp[[2,4,3,l]] = -Sqrt [Expand [temp [ [2 ,4,3 , 1] ] ~ 2] //. ruleqX]; 
temp [[4, 2, 1,3]] = Sqrt [Expand [temp [ [4,2, 1 ,3] ] ~2] //. ruleqX]; 
Expand [temp //. Sqrt [X] -> Y //. X -> Y"2] //. ruleqatop 

]; 

ungradedsigma := 
Table [ 

(- 1) " (Grading [ket [ j ] ] * (Grading [ket [k] ] +Grading [ket [1] ] ) ) * 

sigma[[i, j ,k,l]] , 
{i,M},{j,M},{k,M},{l,M} 

]; 

matrixof ungradedsigma := 
Table [ 

ungradedsigma [ [Ceiling [m/M] , Ceiling [n/M] , 

m+M(l-Ceiling[m/M]), n+M(l-Ceiling[n/M] )]] , 
{m,M-2},{n,M-2} 

]; 

matrixof ungradedsigmainv := 

Simplify [Inverse [matrixof ungradedsigma] ] ; 

ungradedsigmainv := ungradedsigmainv = 
Table [ 

matrixof ungradedsigmainv [ [M(i-l)+k,M(j-l)+l] ] , 
{i,M>,{j,M},{k,M},{l,M} 

]; 

P := P = 

Table [ 

(-1) "Grading [ket [j] ] *KroneckerDelta[i ,1] *KroneckerDelta[k, j] , 
{i,M},{j,M},{k,M},{l,M} 

]; 

Is := Limit [sigma //. Y->Sqrt [-l-q"2+p~ (-2)+q"2 p"2] //. p->q"a, q->l] ; 
(* *) 



^ r\r\ 



CheckRu := CheckRu = 
Module [ 
{ 

temp = 
( 

- (q'u-q'(2a))/(l-q~(u+2a)) Projector [1] 

- Projector [2] 

- (l-q"(u+2a+2))/(q~u-q~(2a+2)) Projector [3] 

) 

}, 

temp = temp //. Qbracket[q_, x_] :> (q"x - q~(-x))/(q - q~(-l)); 
temp = Factor [temp] ; 

temp[[2,l,3,4]] = Sqrt [Simplify [ExpandCtemp [ [2, 1 ,3,4] ] "2] ]] ; 
temp[[l,2,4,3]] = - temp [ [2 , 1 , 3,4] ] ; 

temp [[3, 1,2, 4]] = Sqrt [Simplify [Expand [temp[[3, 1,2,4]] ~2]]] ; 
temp[[l,3,4,2]] = - temp [ [3, 1 , 2,4] ] ; 

temp [[2, 4, 3,1]] = Sqrt [Simplify [Expand [temp [ [2,4,3, 1] ] "2] ]] ; 
temp[[4,2,l,3]] = - temp [ [2 ,4, 3, 1] ] ; 

temp[[3,4,2,l]] = Sqrt [Simplify [Expand[temp [ [3,4, 2, 1] ] "2] ]] ; 
temp[[4,3,l,2]] = - temp [ [3,4, 2, 1] ] ; 

temp = Simplify [temp] ; 

temp = temp //. ((q"(x_) - 1) :> q"(x/2) Q[x/2] (q-q~(-l))); 
temp = temp //. ((q"u-q"x_) :> q"(u/2+x/2) Q [u/2-x/2] (q-q"(-l))); 
temp = temp //. ((q~x_-q~u) :> q~(x/2+u/2) Q[x/2-u/2] (q-q~(-l))); 
temp = temp //. (Q[-a+x_] :> -Q[a-x]); 
Simplify [temp] 

]; 

Format [Q [x_] ] := SequenceForm[" [" , x, "]_q"]; 

(* *) 



1 m 



B.2 LinksGouldlnvariant .m 



(* 



LinksGouldlnvariant . m 

Mathematica code to compute various link polynomial invariants for 
various (oriented) links. We define the positive crossing tensor 
Sigma (here called RMatrix) , its inverse (here called SMatrix) , 
several matrices derived from these, and the 'cup-cap' matrices 
OMMatrix, OPMatrix, UMMatrix, UPMatrix. Requires "Uqgl21.m". 

David De Wit, Louis H Kauffmain, Jon R Links 



04 August 
24 September 
06 October 



1997 — 21 

1997 ~ 04 

1998 ~ 24 



August 1997 
December 1997 
November 1998 



(* DimV is the remge of the indices in the tensors, i.e. the 
dimension of the underlying representation. *) 

DimV [RawBracket] := 2; 

DimV [LinksGould] := 4; 

(* 

UMMatrix [LinksGould] := DiagonalMatrix[{p~2, -p~2, -q~2 p~2, q~2 p~2}] 

UMMatrix [RawBracket] := {{0, I A}, {- I A-(-l), 0>}; 

OMMatrix [Invariant.] := OMMatrix [Invariant] = 
Inverse [UMMatrix [Invariant] ] ; 

OPMatrix [RawBracket] : = OMMatrix [RawBracket] ; 

OPMatrix [LinksGould] := IdentityMatrix [DimV [LinksGould] ] ; 

UPMatrix [Invariant.] := UPMatrix [Invariant] = 
Inverse [OPMatrix [Invariant] ] ; 
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YangBaxterChecker [Invariant _] : = 

Module [ 
{ 

M = DimV [Invariant] , r = RMatrix [Invariant] , temp 

}, 

temp = 
Table [ 
Sum[ 

r[[a,i,b,j]] r[[j,k,c,f]] r[[i,d,k,e]] 

r[[b,j,c,i]] r[[a,d,j,k]] r[[k,e,i,f]] , 
{i,M},{j,M},{k,M> 

], 

{a,M},{b,M},{c,M},{d,M},{e,M},{f ,M} 

]; 

Apply [And, Map[(# == 0)&, Flatten [Expand[temp]]]] 

]; 



RMatrix [RawBracket] := RMatrix [RawBracket] = 
With[ 
{ 

M = DimV [RawBracket] , OM = OMMatrix [RawBracket] , 
UM = OMMatrix [RawBracket] 

}. 

Table [ 

A~(-l) OM[[i,k]] UM[[j,l]] 

+ 

A KroneckerDelta[i, j] KroneckerDelta[k,l] , 

{i,M},{j,M},{k,M},{l,M} 

] 

]; 

RMatrix [LinksGould] := RMatrix [LinksGould] = 

Expand [ungradedsigma /. Y -> Sqrt[p~(-2) - q"2 + p~2 q~2 - 1]]; 

SMatrix [Invaricint_] := SMatrix [Invarieint] = 
Module [ 
{ 

M = DimV [Invarieint] , matrixofr, matrixofs, 
r = RMatrix [Invariant] 

>, 

matrixofr = 
Table [ 

r [ [Ceiling [m/M] , Ceiling [n/M] , 

m+M(l-Ceiling[m/M]) , n+M(l-Ceiling[n/M] )]] , 
{m,M"2},{n,M"2} 

]; 

matrixofs = Inverse [matrixofr] ; 
Table [ 

Expand[matrixof s [[M(i-l)+k,M(j-l)+l]]] , 
{i,M},{j,M},{k,M>,{l,M} 

] 

]; 
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XTensor [Invariaiit_ , X_] : = 
(XTensor [Invariant , X] = 
Switch [ 
X, 

R, RMatrix [Invaricint] , 

RD, XDTensor [Invariant , R] , 

] 

) /; MemberQ[{R, S, RD, SD}, X] ; 



S , SMatrix [Invariant] , 
SD, XDTensor [Invariant , S] 



XLTensor [InvaricLn.t_, X_] : = 
(XLTensor [Invariant, X] = 
With[ 
{ 

M = DimV [Invariant] , Xint = XTensor [Invariant , X] , 

OM = OMMatrix [Invariant] , UM = UMMatrix[Invaricint] 

>, 

Table [ 

Expand [Sum [Xint [[e,d, a, h]] OM[[b,e]] UM[[h,c]], {e,M},{h,M}] ] 
-[a,M>,-[b,M>,-[c,M},-[d,M} 

] 

] 

) /; Meinberq[{R, S>, X] ; 

XRTensor [Invariant., X_] := 
(XRTensor [Invariant, X] = 
With[ 
{ 

M = DimV [Invariant] , Xint = XTensor [Invariant , X], 

OP = OPMatrix [Invariant] , UP = UPMatrix [Invariant] 

}, 

Table [ 

Expand [Sum [Xint [[c,f,g,b]] UP[[a,f]] OP[[g,d]], {f ,M},-Cg,M}]] 
{a,M>,{b,M>,{c,M},{d,M} 
] 

] 

) /; MemberQ[{R, S}, X] ; 

XDTensor [Invaricint_, X_] : = 
(XDTensor [Invariant, X] = 
With[ 
{ 

M = DimV [Invariant] , Xint = XTensor [Invariant , X], 

OP = OPMatrix [Invariant] , UP = UPMatrix [Invariant] 

}, 

Table [ 

Expand [Sum [ 

Xint[[e,f ,g,h]] UP[[a,h]] OP[[g,b]] UP[[c,f]] OP[[e,d]], 
{e,M},{f ,M},{g,M},{h,M} 
]], 

{a,M},-Cb,M},-Cc,M},-Cd,M} 

] 

] 

) /; MemberQ[-[R, S}, X] ; 



Powerof XTensor [Invariant. , power_Integer?Positive , X_] := 
(PowerofXTensor [Invariant , power, X] = 
Module [ 
{ 

M = DimV [Invariant] , Xpminus , Xint = XTensor [Invariant , X] 

}. 
If [ 

power == 1, 

Xint, 

Xpminus = PowerofXTensor [Inveiriant , power - 1 , X] ; 
Table [ 

Expand[Sum[Xint [[a,e,c,f]] Xpminus [ [e ,b,f ,d] ] , {e ,M} , {f ,M}] ] 
{a,M},{b,M},{c,M>,{d,M} 

] 

] 

] 

) /; MemberQ[-[R, RD, S, SD}, X]; 



XUXDTensor [Invariant_, X_] := 
(XUXDTensor [Invariant , X] = 
With[ 
{ 

M = DimV [Invariant] , Xint = XTensor [Invariaint , X] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] , 
XD = XDTensor [Invariant , X] 

>, 

Table [ 

Expand [Sum [ 

Xint[[a,b,e,f]] XD[[g,h,c,d]] OP[[e,g]] UM[[f,h]], 
{e,M},{f ,M},{g,M},{h,M} 
]], 

{a,M},{b,M},-[c,M},{d,M} 

] 

] 

) /; Meinberq[{R, S>, X] ; 

XDXUTensor [Invariant_, X_] := 
(XDXUTensor [Invariant , X] = 
With[ 
{ 

M = DimV [Invariant] , Xint = XTensor [Invariant , X] , 

OM = OMMatrix [Invariant] , UP = UPMatrix [Invariant] , 
XD = XDTensor [Invariant , X] 

>, 

Table [ 

Expand [Sum [ 

XD[[a,b,e,f]] Xint[[g,h,c,d]] OM[[e,g]] UP[[f,h]], 
{e,M},{f ,M},{g,M},-[h,M} 
]], 

-[a,M>,-[b,M>,-[c,M},-[d,M} 

] 

] 

) /; Memberq[{R, S>, X] ; 
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XLXRTensor [Invariant_ , X_] := 
(XLXRTensor [Invariant , X] = 
With[ 
{ 

M = DimV [Invariant] , 

XL = XLTensor [Invariant , X] , XR = XRTensor [Invariant , X] 

}, 

Table [ 

Expand[Suiii[XL[[a,e,c,f]] XR[[e,b,f ,d]] , {e ,M} , {f ,M}] ] , 
{a,M>,{b,M>,{c,M},{d,M} 

] 

] 

) /; Meinberq[{R, S>, X] ; 

XRXLTensor [Invariant_, X_] := 
(XRXLTensor [Invariant , X] = 
With[ 
{ 

M = DimV [Invariant] , 

XR = XRTensor [Invariant , X], XL = XLTensor [Invariant, X] 

}, 

Table [ 

Expand[Sum[XR[[a,e,c,f]] XL[[e,b,f ,d]] , {e ,M} , {f ,M}] ] , 
{a,M>,{b,M>,{c,M},{d,M} 

] 

] 

) /; MeinberQ[{R, S}, X] ; 



XUDUTensor [Invariant., 1 , X_] : = 
(XUDUTensor [Invariant , 1 , X] = 

XTensor [Invariant, X] 
) /; Meinberq[{R, S}, X] ; 

XUDUTensor [Invariant_, N_, X_] := 
(XUDUTensor [Invariant , M, X] = 
With[ 
{ 

M = DimV [Invariant] , XDXU = XDXUTensor [Invariant , X], 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] , 
XUDUTensorLess = XUDUTensor [Invariant , N-2, X] 

>, 

Table [ 

Expand [Sum [ 

XUDUTensorLess[[a,b,e,f]] XDXU[[g,h]] OP[[e,g]] UM[[f,h]], 
{e,M},{f ,M},{g,M>,{h,M} 
]], 

{a,M},{b,M} 

] 

] 

) /; (OddQ[N] && (M >= 3) kk Memberq[{R, S}, X]); 
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AbstractTensor [Invariant., Unknot, x_, y_] := KroneckerDelta[x,y] ; 

AbstractTensor [Invariant., HopfLink, x_, y_] := 
With[ 
{ 

M = DimV [Invaricint] , 

RSquared = PowerofXTensor [Invariant , 2, R] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] 

}, 

Sum [RSquared [[y,x, a, b]] OP[[a,c]] UM[[b,c]], {a,M},{b,M},{c,M}] 

]; 

AbstractTensor [Invariant., Trefoil, x_, y_] := 

With[ 
{ 

M = DimV [Invariant] , 

RCubed = PowerofXTensor [Invariant , 3, R] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invaricint] 

>, 

Sum [RCubed [[y,x,c,d]] OP[[c,f]] UM[[d,f]], {c,M},{d,M},{f ,M}] 



AbstractTensor [Invariant., FigureEight, x_, y_] := 
With[ 
{ 

M = DimV [Invariant] , 

RLRR = XLXRTensor [Invariant , R] , SUSD = XUXDTensor [Invariant , S] , 
UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] 

}, 

Sum[ 

RLRR[[y,a,b,c]] SUSD [ [a, x, c ,d] ] OM[[b,e]] UP[[d,e]], 
{a,M},{b,M},{c,M},{d,M},{e,M} 

] 

]; 

AbstractTensor [Invariant., Cinquefoil, x_, y_] := 
With[ 
{ 

M = DimV [Invaricint] , 

RFif th = PowerofXTensor [Invariant , 5 , R] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] 

}, 

Sum [RFif th[[y,x,c,d]] OP[[c,f]] UM[[d,f]], {c,M},{d,M},{f ,M}] 

]; 

AbstractTensor [Invariant., Septfoil, x_, y_] := 
With[ 
{ 

M = DimV [Invariant] , 

RSeventh = PowerofXTensor [Invaricint , 7, R] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invaricint] 

>, 

Sum [RSeventh [[y,x,c,d]] OP[[c,f]] UM[[d,f]], {c ,M> , {d,M} , {f ,M>] 



AbstractTensor [Invciriajit_, FiveTwo, x_, y_] : = 

With[ 
{ 

M = DimV [Invariant] , ZThree = XUDUTensor [Invariant , 3, R] , 

RSquared = Powerof XTensor [Invarieint , 2, R] , 

UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] 

>, 

Sum[ 

ZThree[[b,c,d,x]] RSquared[[a,b,y,d]] OM[[e,a]] UP[[e,c]], 
{a,M},{b,M>,{c,M},{d,M},{e,M} 

] 

]; 

AbstractTensor [Invariant., WhiteheadLink, x_, y_] := 
Module [ 
{ 

M = DimV [Invariant] , RRRL = XRXLTensor [Invariant , R] , W, 

RDSquared = Powerof XTensor [Invarieint , 2, RD] , 

SSquared = Powerof XTensor [Invariant , 2, S] , 

OP = OPMatrix [Invariant] , UP = UPMatrix [Invariant] , 

UM = UMMatrix [Invariant] 

}, 
W = 
Table [ 

Expeind [Sum [ 

SSquared [ [c, j ,e,f]] RDSquared [ [g , h , i , d] ] OP[[e,g]] UM[[f,h]], 
{e,M>,{f ,M},{g,M},{h,M> 
]], 

{c,M},{j,M},{i,M},{d,M} 

]; 

Sum[ 

RRRL[[a,i,y,b]] W[[c,x,i,d]] OP[[c,a]] UP[[d,b]], 
{a,M},{b,M},{c,M},{d,M},{i,M} 

] 

]; 

AbstractTensor [Invciriant_, SixOne, x_, y_] : = 
Module [ 
{ 

M = DimV [Invariant] , ZThree = XUDUTensor [Invariant , 3, S] , 

SD = XDTensor [Invariant , S] , RRRL = XRXLTensor [Invariant , R] , 
UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] , 
UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] , SOA 

}, 

SOA = 

Table [ 
Sum[ 

SD[[b,c,f ,h]] ZThree[[g,i,d,x]] OM[[f,g]] UP[[h,i]], 
{f ,M>,{g,M>,{h,M},{i,M} 

], 

-[b,M>,{c,M>,-[d,M} 

]; 

Sum[ 

SOA[[b,c,d]] RRRL[[a,b,y,d]] OP[[e,a]] UM[[e,c]], 
{a,M},{b,M},{c,M},{d,M},{e,M} 

] 

]; 
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AbstractTensor [Invciriajit_, SixTwo, x_, y_] : = 

Module [ 
{ 

M = DimV [Invariant] , SRSL = XRXLTensor [Invariant , S] , 

RCubed = PowerofXTensor [Invariant , 3, R] , 
RD = XDTensor [Invariant , R] , 

UM = UMMatrix [Invariant] , DP = DPMatrix [Invariant] , 

UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] , STA 

}, 

STA = 
Table [ 
Sum[ 

RCubed[[e,f ,y,g]] RD[[a,b,c,d]] OM[[c,e]] UP[[d,f]], 
{c,M},{d,M},{e,M},{f ,M} 

], 

{a,M},{b,M},-[g,M} 

]; 

Suin[ 

STA[[a,b,g]] SRSL[[b,h,g,x]] OP[[i,a]] UM[[i,h]], 
{a,M},{b,M},{g,M},{h,M},{i,M} 

] 

]; 

AbstractTensor [Invariant., SixThree, x_, y_] := 
Module [ 
{ 

M = DimV [Invariant] , 

Rint = XTensor [Invariant , R] , Sint = XTensor [Invariant, S] , 
RSquared = PowerofXTensor [Inveiricint , 2 , R] , 
SSquared = PowerofXTensor [Invariant , 2, S] , 

UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] , STA, STB, ST 

}, 

STA = 

Table [ 

Suin[SSquared[[a,b,e,f]] Rint [ [d,e ,y , i] ] , {e,M}] , 
{a,M},{b,M},{d,M},-[f ,M},-[i,M> 

]; 

STB = 
Table [ 

Sum [Sint [[b,c,g,h]] RSquared [ [f ,g, i ,x] ] , {g,M}] , 
{b,M},{c,M},{f ,M},{h,M},{i,M} 

]; 

ST = 
Table [ 

Sum[STA[[a,b,d,f ,i]] STB[[b,c,f ,h,i]] , {b,M>,-[f ,M>,-[i,M}] , 
-[a,M},-[c,M},-[d,M},{h,M} 

]; 

Sum[ 

ST[[a,c,d,h]] DM[[j,d]] UP[[j,h]] OM[[k,a]] UP[[k,c]], 
{a,M},{c,M},-[d,M>,-[h,M},{j,M},-Ck,M} 

] 

]; 
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AbstractTensor [Invciriaiit_, SevenTwo, x_, y_] : = 

With[ 
{ 

M = DimV [Invariant] , ZFive = XUDUTensor [Invariant, 5, R] , 

RSquared = Powerof XTensor [Invarieint , 2, R] , 

UP = UPMatrix [Invariant] , OM = OMMatrix [Invariant] 

>, 

Suin[ 

ZFive[[b,c,d,x]] RSquared[[a,b,y,d]] OM[[e,a]] UP[[e,c]], 
{a,M},{b,M},{c,M},{d,M},{e,M} 

] 

]; 

AbstractTensor [Invar iaiit_, EightSeventeen, x_, y_] := 
Module [ 
{ 

M = DimV [Invariant] , EA, EB, EC, ED, ES, 

Rint = XTensor [Invariant , R] , Sint = XTensor [Invarieint, S] , 

RSquared = Powerof XTensor [Invariant , 2, R] , 

SSquared = Powerof XTensor [Invariant , 2, S] , 

UM = UMMatrix [Invariant] , OP = OPMatrix [Invariant] 

>, 

EA = 
Table [ 

Expand [Sum [RSquared[ [a, b,g,d]] SSquared[[c,g,e,f]] , {g,M}]], 
{a,M},{b,M},{c,M},{d,M},-Ce,M},{f ,M} 

]; 

EC = 
Table [ 

Expand [Sum [Sint [[d,k,f,l]] Rint [ [b,m,k,n] ] , {k,M>]], 
{b,M},-[m,M},-[d,M},{n,M},{f ,M},{1,M} 

]; 

ED = 
Table [ 

Expand [Sum [Sint [[n,o,l,j]] Rint [ [m,h,o, i] ] , {o,M}]], 
-[m,M},{h,M},{n,M},-[i,M},{l,M},{j,M} 

]; 

EB = 
Table [ 

Expand [Sum [ 

EC[[b,m,d,n,f ,1]] ED[[m,h,n,i,l, j]] , 
{l,M},{m,M},{n,M> 
]], 

{b,M>,{h,M>,{d,M},{i,M},-[f ,M>,{j,M} 

]; 

ES = 
Table [ 

Expand [Sum [ 

EA[[a,b,c,d,e,f]] EB[[b,h,d,i,f , j]] , 
{b,M},{d,M},{f ,M> 
]], 

-[a,M>,-[h,M>,-[c,M},-[i,M},{e,M},{j,M} 

]; 

Sum[ 

ES[[y,x,c,i,e,j]] OP[[c,r]] UM[[i,r]] OP[[e,s]] UM[[j,s]], 
{c,M},{e,M},{i,M},{j,M},{s,M>,-Cr,M} 

] 

]; 
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AbstractTensor [Invciriajit_, NineFortyTwo , x_, y_] : = 

Module [ 
{ 

M = DimV [Invariant] , SCubed = Powerof XTensor [Invarieint, 3, S] , 
RDSquared = Powerof XTensor [Invarieint , 2, RD] , N, 

KURD = XUXDTensor [Invariant , R] , SDSU = XDXUTensor [Invariant , S] , 
OM = OMMatrix [Invariant] , UP = UPMatrix [Invariant] , 

OP = OPMatrix [Invariant] 

>, 
N = 

Table [ 

Expand [Sum [ 

RDSquared[[a,b,c,d]] SCubed[[e,f ,g,h]] OM[[c,e]] UP[[d,f]], 
{c,M},{d,M},{e,M},{f ,M} 
]], 

{a,M},{b,M},-[g,M},{h,M} 

]; 

Suin[ 

N[[a,b,y,h]] SDSU[[b,i,h,j]] RURD[[k,x,i,m]] UP[[in,j]] OP[[k,a]], 
{a,M},{b,M},{h,M},{i,M},{j,M},{k,M},{in,M} 

] 

]; 

AbstractTensor [Invariant., TenFortyEight , x_, y_] := 
Module [ 
{ 

M = DimV [Invariant] , TA, TB, TT, 

Rint = XTensor [Invariant , R] , 

SSquared = Powerof XTensor [Inveiriajit , 2 , S] , 

SCubed = Powerof XTensor [Invariant , 3 , S] , 

RFourth = Powerof XTensor [Invariant , 4, R] , 

UP = UPMatrix [Invariant] , OP = OPMatrix [Invariant] , 

OM = OMMatrix [Invariant] , UM = UMMatrix [Invariant] 

>, 

TA = 
Table [ 

Expand [Sum [SSquared [ [a , b , y , f ] ] RFourth [ [f , d , g , h] ] , {f , M}] ] , 
{a,M},{b,M},{d,M},{g,M},-Ch,M} 

]; 

TB = 
Table [ 

Expand [Sum [SCubed [[b,c,d,e]] Rint [ [e,x,h, i] ] , {e,M}]], 
-[b,M>,-[c,M>,-[d,M},-[h,M},{i,M} 

]; 

TT = 
Table [ 

Expand [Sum [ 

TA[[a,b,d,g,h]] TB[[b,c,d,h,i]] , 
{b,M},{d,M>,{h,M} 
]], 

-[a,M>,{c,M>,-[g,M},-[i,M} 

]; 

Sum[ 

TT[[a,c,g,i]] OM[[j,a]] UP[[j,c]] OP[[g,k]] UM[[i,k]], 
{a,M>,{c,M},{g,M},{i,M},{j,M>,{k,M} 

] 

]; 
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AbstractTensor [Invciriajit_, Pretzel [pp_, qq_, rr_] , x_, y_] : = 

Module [ 
{ 

M = DimV [Invariant] , PI, P2, P3, 

Zp = XUDUTensor [Invaricint,pp,R] , Zq = XUDUTensor [Invaricint,qq,R] , 
Zr = XUDUTensor [Invariant, rr,R] , 

OM = OMMatrix [Invariant] , UP = UPMatrix [Invariant] 

}, 

PI = 

Expand [Table [ 

Sum[Zp[[a,b,y,e]] * Zq[[b,c,e,f]] , {b,M},{e,M}] , 
{a,M>,{c,M},{f ,M} 

]]; 

P2 = 

Expand [Table [ 

Sum[Pl[[a,c,f]] * Zr[[c,d,f ,x]] , {c,M},{f ,M}] , 
{a,M},{d,M> 

]]; 

P3 = 

Expeind [Table [ 

Sum[P2[[a,d]] * OM[[g,a]], {a,M}] , 
-[d,M>,-[g,M> 

]]; 

Expand[Sum[P3[[d,g]] * UP[[g,d]], {d,M} , {g,M}] ] 
] / ; ( (OddQ [pp] && OddQ [qq] && OddQ [rr] ) && (3 <= pp < qq < rr) ) ; 

KTATensor [Invarieint_] : = 
Module [ 
{ 

M = DimV [Invariant] , SD = XDTensor [Invariant , S] , 

SSquared = PowerofXTensor [Invariant , 2, S] , 
OP = OPMatrix [Invariant] , UM = UMMatrix [Invariant] , 

RURD = XUXDTensor [Invariant , R] , temp 

}, 

temp = 
Table [ 

Expand [Sum [ 

SSquared [[d,s,f,g]] SD [ [h, i , e , c] ] OP[[f,h]] UM[[g,i]], 
{f ,M},{g,M},{h,M>,{i,M} 
]], 

{d,M},{s,M},{e,M},{c,M} 

]; 

Table [ 

Expand [Sum [RURD [ [a, d,b,e]] temp [ [d, s , e , c] ] , ■[d,M},-Ce,M}] ] , 
{a,M},{s,M},{b,M>,{c,M} 

] 

]; 
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KTAITensor [Invariant _] : = 
Module [ 
{ 

M = DimV [Invariant] , SD = XDTensor [Invariant , S] , 

SSquared = Powerof XTensor [Inveiricint , 2, S] , 
OP = OPMatrix [Invariant] , UM = UMMatrix [Invariant] , 

KURD = XUXDTensor [Invariant , R] , temp 

}, 

temp = 
Table [ 

Expand [Sum [ 

SSquared[[a,d,f ,g]] SD[[h,i,b,e]] OP[[f,h]] UM[[g,i]], 
{f ,M},{g,M},{h,M},{i,M} 
]], 

{a,M},{d,M},{b,M>,{e,M} 

]; 

Table [ 

Expand [Sum [temp [ [a, d,b,e]] KURD [ [d, s , e , c] ] , ■Cd,M},{e,M}] ] , 
{a,M},{s,M},{b,M},{c,M} 

] 

]; 

KTCTensor [Invariant _] : = 
Module [ 
{ 

M = DimV [Invariant] , KTB, 

OP = OPMatrix [Invariant] , OM = OMMatrix [Invariant] , 

UP = UPMatrix [Invariant] , UM = UMMatrix [Invariant] , 

SLSR = XLXRTensor [Invariant , S] , Rint = XTensor [Invariant, R] , 
RDSquared = Powerof XTensor [Invariant , 2, RD] , 
SUSD = XUXDTensor [Invariant , S] 

}, 

KTB = 
Table [ 

Expand [Sum [ 

Rint[[d,a,b,c]] RDSquared [[l,m,f,n]] SUSD[[a,e,n,g]] * 
OP[[b,l]] UM[[c,m]] , 
{a,M},{b,M},{c,M},{l,M},{m,M},-Cn,M} 
]], 

-[d,M>,-[e,M>,-[f ,M},-[g,M} 

]; 
Table [ 

Expand [Sum [ 

KTB[[d,e,f ,g]] SLSR[[h,i,j,k]] OM[[f,h]] UP[[g,i]], 
{f ,M>,{g,M>,{h,M},{i,M} 
]], 

{d,M},{e,M},{j,M},{k,M} 

] 

]; 
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AbstractTensor [Invciriaiit_, KT, x_, y_] : = 

With[ 
{ 

M = DimV [Invariant] , KTC = KTCTensor [Invariant] , 

KTA = KTATensor [Invariant] , UP = UPMatrix [Invariant] , 
OP = OPMatrix [Invariant] , OM = OMMatrix [Invariant] 

>, 

Suin[ 

KTA[[a,x,b,c]] KTC[[d,e,j,k]] * 
OM[[b,d]] UP[[c,e]] OP[[a,j]] UP[[k,y]], 
{a,M},{b,M},{c,M},{d,M},{e,M},-Cj,M},{k,M} 

] 

]; 

AbstractTensor [Invariant., KTI, x_, y_] := 
With[ 
{ 

M = DimV [Invariant] , KTC = KTCTensor [Invariant] , 

KTAI = KTAITensor [Invariant] , UP = UPMatrix [Invariant] , 
OP = OPMatrix [Invariant] , OM = OMMatrix [Invariant] 

}, 

Suin[ 

KTAI[[a,x,b,c]] KTC[[d,e, j ,k]] * 
OM[[b,d]] UP[[c,e]] OP[[a,j]] UP[[k,y]], 
{a,M},{b,M>,{c,M},{d,M},{e,M},{j,M},{k,M} 

] 

]; 

(* *) 
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KnotList = 
{ 

Unknot, HopfLink, Trefoil, FigureEight , 
Cinquefoil, FiveTwo, WhiteheadLink, 
SixOne, SixTwo, SixThree, Septfoil, SevenTwo, 
EightSeventeen, NineFortyTwo , TenFortyEight , KT, KTI 

}; 

LinkPolynomial :: usage = "LinkPolynomial [Invariant , Link] evaluates the 
link polynomial invariant of narnie Invariant for the knot described by 
Link. Invaxiant must be one of ' RawBracket ' , or 'LinksGould' , and Link 
must be one of the elements of 'KnotList' or 'Pretzel[pp, qq, rr] ' . " 

Attributes [LinkPolynomial] = {Listable} ; 

LinkPolynomial [RawBracket , Link_] := 
LinkPolynomial [RawBracket , Link] = 
With[ 
{ 

M = DimV [RawBracket] , 

OM = OMMatrix [RawBracket] , UP = UPMatrix [RawBracket] 

>, 

Print ["Building LinkPolynomial [RawBracket, ", Link, "]."]; 
Expand [Sum [ 

AbstractTensor [RawBracket , Link, x, y] OM[[z,y]] UP[[z,x]], 
{x,M},{y,M},{z,M} 

]] 

]; 

LinkPolynomial [LinksGould, Link_] := 
LinkPolynomial [LinksGould, Link] = 
( 

Print ["Building LinkPolynomial [LinksGould, ", Link, "]."]; 
Expand [AbstractTensor [LinksGould, Link, 1, 1]] 

); 

(* *) 

Attributes [LGCcintDetectMoninvertibilityQ] = {Listable}; 

LGCcintDetectNoninvertibilityCl [Link_] : = 
With[ 

{poly = LinkPolynomial [LinksGould, Link] /. p -> q"a}. 
Expand [poly - (poly /. a -> -(1+a))] === 

]; 

Attributes [LGDetectsChiralityQ] = {Listable}; 

LGDetectsChiralityQ [Link_] : = 
With[ 

{poly = LinkPolynomial [LinksGould, Link] /. p -> q~a}, 
PowerExpand [poly - (poly /. q-> q"(-l))] =!= 

]; 

(* *) 
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Listof Pretzels [N_] := 
Flatten [ 

Table [Pretzel [pp.qq.rr] , {pp,3,N,2}, {qq,pp+2,N,2>, {rr,qq+2,N,2>] 

]; 



LGCantDetectNoninvertibilityof PretzelsQ [N_] : = 
LGCantDetectNoninvertibilityQ [Listof Pretzels [N] ] ; 

LGDetectsChiralityof PretzelsQ [N_] := 
LGDetectsChiralityQ [Listof Pretzels [N] ] ; 



(* 



*) 



Writhe [Pretzel [pp_ , qq_, rr_]] := pp + qq + rr; 



Writhe [Unknot] = 

Writhe [Hopf Link] = 2 

Writhe [FigureEight] = 

Writhe [FiveTwo] = 5 

Writhe [SixOne] = -2 

Writhe [SixThree] = 

Writhe [SevenTwo] = 7 
Writhe [NineFortyTwo] = -1 

Writhe [KT] = -2 



Writhe [Trefoil] = 3 

Writhe [Cinquef oil] = 5 

Writhe [WhiteheadLink] = 2 

Writhe [SixTwo] = 2 

Writhe [Septf oil] = 7 
Writhe [EightSeventeen] = 

Writhe [TenFortyEight] = 

Writhe [KTI] = -2 



Attributes [JonesPolynomial] = {Listable}; 

JonesPolynomial [Link_] : = 
With[ 
{ 

bu = LinkPolynomial [RawBracket , Unknot], 
b = LinkPolynomial [RawBracket , Link] 

}, 

Expand [(-A) "(-3 Writhe [Link] ) Expand [Factor [b/bu]]] /. A -> t~(-l/4) 

]; 

Tableof JonesCoef f s [Listof Links_] : = 
Module [ 
{ 

J = JonesPolynomial [Listof Links] , powers 

}, 

powers = -Exponent[J, t"(-l)]; 
Transpose [ 
{ 

Listof Links , 

powers , 

CoefficientList [Expand[t~ (-powers)] J, t] 

} 

] 

]; 



(* 



*) 



lie 
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